Problèmes aux limites de Laplace - secteur de disque en coordonnées polaires — p 
Etude des problèmes aux limites de l'équation de Laplace sur des secteurs de disque. 
Ce sont des problèmes de deux types que nous allons étudier 


Le premier est homogène dans la dimension angulaire 
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On rappelle que cela peut poser des problèmes si l'on souhaite prendre en compte des conditions 
mixtes, car en théorie elles s'écrivent sous cette forme : 
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et ce n'est qu'une aa raoi que de considérer les termes comme fixes. Dans tous les cas en 
donnant les valeurs limites de ces coefficients des conditions aux limites, on revient aux cas 
homogènes de Dirichlet et Neumann, qui ne posent pas de problème théorique majeur. 


Le second est homogène dans la dimension radiale : 
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Pour le problème homogène radiale, il n'y a pas de souci pour la formulation des conditions 
homogène mixte. 
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Pour ces deux types de problème, la question demeure la suivante, le système séparé et homogène 
dans la dimension choisie permet-il d'exhiber une base de fonctions propres et toutes les bonnes 
propriétés d'un système de Sturm-Liouville et ainsi d'envisager sereinement le développement en 
série des solutions. Nous verrons que c'est même le cas pour la dimension radiale, dans le cas d'un 
secteur annulaire. 
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Exemple: Secteur de disque plein (r,0) soumis à des conditions aux limites de Dirichlet 
homogènes en 0 


Occupons maintenant d'un secteur de disque plein soumis à des conditions de Dirichlet sur 
chacune de ces frontières. Sans restreindre le problème on peut fixer le rayon du disque à 1. Pour 
des raisons de plus grande simplicité de la solution cherchée on place le secteur de disque posé sur 
l'axe des abscisses. Toute position différente trouverait sa solution par un simple changement de 
variable en @ correspondant à une rotation adéquate. 
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On rappelle que les solutions de l'équation séparée en (r,0) sont de la forme 


T(r,0)=e+ fin(r)+ X (e,r” +d,r * (a, Cos(2,0)+ b,Sin(2,0)) 


n=1,00 


Ici il faut bien faire attention à revenir à un spectre général de valeurs propres selon la variable 0 
au delà des valeurs entières comme dans les cas précédent qui ne sont liées qu'à une géométrie de 
disque entier, et non sectoriel. 


La condition de nullité en 6=0 donne an=0, la condition de finitude de la solution donne dn=0 et 
f=0. De la configuration des conditions aux limites en r=0, réunion de deux conditions nulle sur les 
tranches du secteur, on déduit facilement que la valeur de la solution en r=0 est nulle également. 
Donc e=0. La solution prend la forme plus simple : 


T(r,0)= $ b,r™ Sin(2,0) 
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Appliquons maintenant la conditions aux limites en O=a 
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Si les conditions aux limites porte sur l'intervalle [0,02] alors recherchons la solution sous la forme 
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Il suffit finalement d'appliquer la condition aux limites en r=1, pour trouver la solution du 
problème, en utilisant les relations d'orthogonalité des fonctions propres @(8,. ); Sur l'intervalle 
0e[0,,8,]- [0,27] soit: 
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ce qui donne : 
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La détermination des valeurs propres correspond au cas où l'angle a est aigu. Dans ce cas la 
condition de régularité du gradient est assurée. Lorsque l'angle est supérieur à n (angle obtu), la 
finitude du gradient n'est pas assurée, mais par un raisonnement sur la solution construite et son 
extension au domaine complémentaire dans le disque, soit le secteur de l'angle complémentaire 
aigu B=2 n- a, on peut montrer que la solution trouvée possède bien la régularité du gradient 
requise, car elle peut également en théorie être construite à l'aide des fonctions et valeurs propres 
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du problème complémentaire P . Le but n'est pas d'exhiber cette 

construction mais de vérifier la propriété de régularité requise. 
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On ne peut étendre la technique de séparation des variables lorsque les conditions aux limites sont 
inhomogènes en 0 avec des conditions très générales angulaires, sur un secteur de disque plein. 
Hélas, il n'existe pas de système de fonctions propres radiales issues de la séparation de Lamé qui 
possèdent les propriétés requises : 
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L'illustration de tels cas montre bien que l'utilisation de la séparation des variables même quand 
elle est possible ne conduit pas toujours à la solution recherchée. Précisons pourquoi ces fonctions 
propres n'existent pas : on rappelle que la séparation des variables de l'équation de Laplace en 
coordonnées polaire, donne : 
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La constante peut être soit positive, soit négative, soit nulle et les valeurs doivent être choisit en 
cohérence avec la configuration géométrique. Dans le cas x , les solutions sont : 

r (r R'(r))-ÆR(r)=0= R(r) = ar” +br” 

©" (0) +200) = 0 = @(0)= e Cos(20)+ f Sin(20) 

Dans le cas + les solutions sont 
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Et dans le cas d'un disque complet on doit écarter le choix -7 par le fait que les solutions 


angulaires n'étaient pas périodiques. De même pour une disque plein, ou toute autre géométrie 
qui inclut dans le domaine l'origine, la solution de la forme 


e ih —i A . DN . . . . . . . 
R(r)=ar" +br" était à proscrire impliquant une oscillation rapide de la solution autour de 
l'origine. Dans ces mêmes configurations, il n'est donc pas possible d'exhiber une classe de 
fonctions propres notamment dans la dimension radiale. 
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Dans ce cas il faut alors se tourner vers des techniques de calcul plus complexe, par exemple dans 
la catégorie soit des transformations intégrales et du calcul opérationnel, soit les transformations 
conformes du plan complexe. Peu de cas de l'utilisation de ces techniques sont rapportés dans la 
littérature scientifique, on peut noter l'utilisation des transformations de Meilin dans le plan 
complexe, ainsi que les transformations Sinus et Cosinus logarithmique. La technique consiste, 
sans entrer dans le détail de complexité, à utiliser une transformation intégrale vers le plan 
complexe d'une nouvelle variable afin de "simplifier" le calcul de la solution transformée. Puis de 
revenir au calcul de la solution dans le plan initial par la transformée inverse. Cette dernière étape 
en en général plus complexe, car elle met en oeuvre des techniques d'intégrales de contour dans le 
plan complexe par la méthode de Cauchy. 


Mais attention n'allons pas trop vite en besogne, il y a un cas où l'on peut construire une base de 
fonctions propres dans la dimension radiale, celui où le domaine n'inclut pas l'origine, et comme 
dans la dimension angulaire, il s'agit d'un secteur d'angle, il n'y a alors plus la nécessité de 
respecter la périodicité en Ÿ ! Mais nous aborderons ce cas plus loin dans ce document 


Quelques cas simples de solutions par séparation des variables aux problème aux limites 
inhomogènes angulaires 


Revenons pour l'instant à un secteur de disque plein...avec les problèmes suivants homogène enr 
et inhomogène en 9, avec une valeur constante en condition limite, il est possible de construire la 
solution à l'aide du principe de superposition : 
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En effet, utilisons la solution triviale du problème aux limites de Dirichlet dans un domaine non 
bornée (secteur d'angle a dans le plan) : 
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condition aux limites du second problème devenu cette fois-ci homogène, mais avec une condition 
aux limites qui dépend de l'angle 0 : 
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Ce qui nous permet d'introduire justement dans cet exemple une extension à des problèmes plus 
généraux : 
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Revenons alors aux problèmes décrit sur un secteur de disque, qui conduit aux sous problèmes 
homogènes suivants : 
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La solution formelle du problème homogène définie plus haut est la suivante : 


| nei y 


a n=1,00 





nī 


An 
(= b Sin(2,0) 


b, -fao f, (0) Sin(A,0) À, = 


ce qui donne : 























Comme 
i Sin(2,0)-À,0 Cos(1,0)f ~% 2-470? ,0)+ 22,0 Sin(À,0 
| d0 0 Sin(2,0) _[sin(,8) a os(2,0)] [d0 0° Sin(2,0) = ( 26° Cost, H 10 Sin(å, f 
0 À; 0 À 
à 3[4,°0? — 2 ]Sin(4,0)+ 4,0470? -6]Cos(2,0 
fao 0 Sin(,0) K n in( n )+ ol n ) os( n E 
0 An 
4 44 0(2 0? -6 lSin(4 0)-(4 "04 -122 70? +24 4,0 
[ d0 6* Sin(2,0) l : b. AR) b, : heu } 
0 À, 
0 i 0 
fo(0)=—-1 fa(0)=-— 
a a 
1$% , Poe 1$ ; 
b, =— | d0 0 Sin(2,0) — | d0 Sin(4,0) b, =-— | d0 0 Sin(2,0) 
a 0 0 a 0 
f d0 Sin(2,0) = CCE) f d0 Sin(2,0) = Le 
0 An 0 An 
f 40 0 Sin(2,0) RAF [d0 0 Sin(2,0) Pa 
0 À 0 À 
paa D _(-6D°)__ 1 ba (-1) 
À, À, À, QM À, 
jee parr 
a a 
0° 2 Ni 8 2 sqa 
da — r r z 
T(r,0)= z Sin(A,0 T(r,0)=—+ Sin(A,0 
(r,0) z z|) 2, in(,0) (r,0) = PHA à, in(,0) 
À À 
a-0 2 r |” Sin(A,0) 0 2 r |" =D" Sin(A,0) 
T ,0 = = i T ,0 = L 
8) a zy|z) n Pa a DAA n 
T(r,0) „fini T(r,0)._. =0 T(r,0) „fini T(r,0). =0 








T(r,0),,=1 T(r,0),_, =0 




















(PED =0 TE 


Problèmes aux limites de Laplace - secteur de disque en coordonnées polaires — p 


On construit donc par principe de superposition la solution du problème suivant avec trois valeurs 
fixes différentes aux limites : 
AT(r,0)=0 T(r,0) fini 
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Construisons par le même procédé, la solution des problèmes : 





AT(r,8)=0 T(r,0)_ fini AT(r,0)=0 T(r,0)_ fini 
ou |T(r,0) = 
T(r,0), =1 Tre =0 T(r,0), =0 T(r,0),. =1 


En effet, utilisons la solution triviale du problème aux limites de Dirichlet dans un domaine non 
bornée (secteur d'angle a dans le plan) : 


T eige (r, 0)= 0 T iedee (r, o) fini Traza (r, 0)= 0 T sedge (r, o) fini 
kegi ee En =1 EE (r,0),. = 0 ou T wedge Oa = = 0 Poe ro) = 
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Ce qui donne les deux solutions : 
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Donnons alors la solution pour les problèmes généraux suivants : 
AT(r,0=0 T(r,0).., fini AT(r,0=0 T(r,0).., fini 
r=l, = fo (0) ou T(r, o), = Ío (8) 
T(r,0),, = Ta TE Opa =0 T(r,0),., =0 T(r,0),, =To 
AT seage t,0)=0 Tuag (rO| _ fini AT vse (T,0)=0 Tragel, O)| Fini 
T sedge T> 0), L eg 0). = ou T sedge T> o) (SN Ta D = To 
8-6, |. F -0 p-r, 28 
=> T edge (f> 0)= T, i 0, ne Ta 0, zg | T sedge F> )= Ty, 6, 6. 
La solution peut s'écrire T(r:0)= Trege (t 0)+ Teeneous (Y0) 


, avec un second problème devenu cette 
fois-ci homogène, avec une condition aux limites modifiée dépend de l'angle 8 : 
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AT 6, -0 | _ nī 
Se er) c, -fao[ n0- É 2) sre. (0-0) dei er) c, ff n0- E 
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gl A, 
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Soit pour le premier problème la solution : 
AT(r,0=0 T(r,0)., fini 


T (0). = /9(0) 
TEA) =T, T(r,0),, = 
À z nī 


i (8,-0,) 


b, = | d0 f,(0) Sin, (0-8) 








E rŸ" Sina (0-0). 2 Pre : 
T(r:0)= 27,5 [2) . Et b,Sin(2,(0 -0,)) 


fa 0, -0, m n=1,0 
Et pour le second problème la solution : 
AT(r,0=0 T(r,0) fini 


2, = Jo(0) 
T(r,0),., = 0 TO) oo, =T; 





0 nT 
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b, = d0 f,(0) Sin(à,(8-0,) 








(0, -0,) T 
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An i À 
T(r,0)=T, 2 — RS B — Sin(4,(8-8,)+,— Y f | b,Sin(à, (0-0) 
2 =1, > 





Problèmes aux limites de Laplace - secteur de disque en coordonnées polaires — p 142 


Finalisons par principe de superposition pour le problème suivant : 
AT(r,0=0 T(r,0)., fini 


TEO), z Jo (0) 
T(r,0),., = To, PO y = Lo 


2 


Lg) = [40 LE) Sn(2,(6-8) 








À 

0,-0 0-0, 2 r |" Sin(A,(8-0,)) 

T(r,0) =T, = +T, L+ i LOC, =T, )+ 
( ) GN 0, -—9, 6, 0, -0, T DIE | n (c ) CA z) 


r: 


TE (= bsna) 


(6, m 0) n=1,00 l 


Exemple: Secteur de disque plein (r,0), illimité dans le plan, soumis à des conditions aux limites 
de Dirichlet. 


Abordons maintenant un problème sur un domaine illimité qui peut également être résolu par 
séparation des variables. Il s'agit du problème suivant : 
AT(r,0=0 T(r,0)., fini 


T OSFSE 
TO Ooo = 0 l <r<o 
T, O<r<l, 
TE Ooo, = 0 l <r<o 


Ramenons ce problème à l'étude des problèmes homogènes : 


Le problème est équivalent à deux sous problème en séparant l'espace en deux sous domaines Q: 
et Q, définit comme suit : 
Q,={0<r<1,0<0<0,} Q, ={,<r<%,0<0<0,} 


T(r,0) =T, (r,0)+To (r,0) 
Ta, 0) = Ta CO) =D 
i, CD N aio (0), = 0 
Continuité To (1,0) = To, (l-0) 


Problèmes aux limites de Laplace - secteur de disque en coordonnées polaires — p 143 


Le problème sur Q, est homogène et la solution s'écrit sous la forme (voir section domaines non 
bornés) : 


CA 
= Cu = [dx f,(6) Sin(4,80) 
0 


0 


2 LA 
To(,0)= = Y e(t) Sin(2,0) 


0 n=1,0 





(2n+1)r 4 É j Sin(2,0) 
0)=1> 4, =“ T(r,0)= r n 
St fol ) n 6, (r ) 2 A 2n+1 
4 5 Sin(,0) _, 
Rue. 2n+1 


Le problème sur Q, peut également être ramené à un problème homogène en séparant en deux 
sous-problèmes, comme suit : 

Q,={0<r<1,0<0<06,} Q,=f, <r<x0<0<0,} 

To, (r,0) = Ta, (r,0)+ To, (r,0) 


Ta, (0), a (0), =h To, CO) =h 


To, (0), = Toa 0), 50 Ta O -+T O | 
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Le premier sous-problème admet la solution triviale : T=T, et le second le développement classique 
sur le secteur de disque borné, soit : 

nī 
À,=— 

CA 


À 
2 Une 
Ta 0)=D+— > af] Sin(2,0) 
0 n=1,0 r 

Il s'agit maintenant de s'assurer des conditions de continuité de la solution, et de continuité de la 
dérivée première radiale à la frontière des deux domaines Q; et Q,, en utilisant la décomposition 
de l'unité par le fonctions propres. Ce faisant on aboutit à une système assez simple d'équations 
linéaires suivant : 











4 Sin(A,0) 2 n\Sin(A,0) nm nm 
nicis] 1—(-1 n Aris 0, = 
T e 2n+1 = 2 l r) n " 6, se. À, 
À, 
2 „\Sin(A,0) 2 = À b ee 
STOR Z h-e) EA, y ARE) sc, 
n=1,00 n=1,00 r 





2 n\, À,b, | Sin(A,0 2 À 
TaD D (nh-e y AE GO 70,02 I LE sinap 


n=1,%0 n=1,00 








2 2 
oTo (l.,0) 2 2 > A, b, Sin(À,0) Lii ie 2 > l À, Ca Sin(A,0) 
y 


ôr m n=1,% n r m n=1,0 n r 


OTG a g 











ôr ôr n Ca 
To, (L.,0) z To (l,,0) > (7 (1-61) à 2.) Ac, €, = nG ) 
n n n 22, 
Ca = T; Ge 1) 90 = -b, 
n 27 
La solution du problème s'écrit donc : 
Re (2n+1)r 
0, 





T no» 2N+1 


US DE > — z | sao) 





2 ME” 
Tat DRS, Le Sin(À,0) 


T now 2 LT 
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Admettons que 8,=n, alors il vient : 
À,=2n+1 


Ta (r,0)= rfi- > B sneno 


T noo 2n+1\1, 








T..(r,0)=T, À > : OI Sin((2n +1)0) 
F. 


T aoo 2N+1 


Sin((2n + 1) = Cos(nx)=(-1)" 


n c : ( 1) 2n+1 
e 2 y CN ; | 


epig H 


noo 2n+1 








T 
Toa r>) y b Si 
Posons gat > 2 aje G 
I Ty T noo 2nN+1 


2n+1 
X 





r 


l+ix 


—] 1 
ED i x”! est le développement de Taylor dela fonction F Lou = 


n=0,00 2n+ Le 





z) et Arctan(x) 
ix 
d'où Profil de T(r,0) en 6 = 


-ży st) n [1-2 Arean) O<x<1 
T 


T oo 2n+1 
roo » xCD LT < Arctan —| 1<x<o 
T oo 2n+1 m X 


1 
Comme pour xe R, Arctan(x)+ Areta +) = - 
x 


>T(r, D5 1,2 Aretan +) - T(r D5 2 Arora £) O<r<o 
r 


Mais pour éviter la singularité enr =0 on préfère écrire 


nh: Aa Arcun t) OSFSE 
Tr, = a | 


2 2 L 
T, — Arctan >| l,<r<Ħæœ 
m r 
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Exemple: Secteur de disque plein (r,0), illimité dans le plan, soumis à des conditions aux limites 
de Dirichlet. 


Continuons sur la lancée du problème précédent, toujours sur un domaine illimité qui peut 
toujours être résolu par séparation des variables. Il s'agit du problème suivant : 
AT(r,0=0 T(r,0)., fini 
0 O<7r 
T(r,0),., = T(r,0), 41, Lier 
0 1, <r<o 
On ramène ce problème à l'étude des trois problèmes homogènes sur trois sous-domaines Q; , Q; 
et Q; définit comme suit : 
Q,={0<r<1,,0<0<0,} Q,={fl,<r<1,0<0<0,} Q,={,, <r<ow,0<0<0,} 
T(r,0)=T, (r,0)+70,(r,0)+To,(r,0) 
Continuité 


To, (0). E To, (r,0)., 





r CA E N GS 


Continuité de la dérivée 




















aTa (rO| Th, (r,0) 

or r=l,; í ôr fel 
aTa, (r,O| OTa(r,6) 

or r=l;z or r=l,2 


Le problème sur Q, est homogène et les solutions sur Q, et Q; s'écrivent sous la forme (voir ci- 
avant) : 

nm 
A, =— 

CA 


T(r,0)=— 2 AÎZ J Sin(À,0) 


One 


T,9= 2 Ÿ D.(° Le) | Sin(,0) 


0 n=1,0 


Le problème sur Q; suit la forme des problèmes aux limites avec deux fonctions limites définies 
sur le rayon SRE Ir1 et extérieur |r2 : 


re B, -fas fI (O) Sin(4,0) C, = jas fE (0) Sin(A,0) 


Ta2(r,0) = To + os (l a) iza (GEGI 
GEGE A 





Sin(å,0)+ 





Sin(A,0) 
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Il s'agit maintenant de s'assurer des conditions de continuité de la solution, et de continuité de la 
dérivée première radiale à la frontière des deux domaines Q; et Q,, en utilisant la décomposition 
de l'unité par le fonctions propres. 


Bore NaS (ESP PE 4 y> Sin(,0) 
Le n 


0 n=1,00 n m n=0,00 2n+1 
Ce faisant on aboutit à une système assez simple des 2 équations linéaires suivantes pour la 
continuité CO: 


4 28, +00) Lo 4,78, -M-CN) = =U 
D, +, D HO) 


Notons 








n=1,0 
































l 


rl 


À 
28,[le | i 
„=B, + - ln — + 2 = er. (fe) +2C',=U 
ARIANE fn 
A) E 
À, 
zca] 2, 
D, =-C, + B; -+ ln & 2B',+2C', fl) =-U 
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Ce qui donne par évidence : 
B', = =C © B, = C, 


__ 0 o an) ga ñt-t) 


n l 4 ’ An n I À 
SEE Sue) AE 
La La La 

















ñ m Àn An 
pega aiey = À, = ni- G) La + Laz —9 
À, I is L La 
22, (2) a 


La solution complète, sachant que la série ne présente que des termes de développement impairs, 
revêt la forme : 
Po (2n+1)r 

n 0, 








À An 
E 
OE 2 3, a RE z 
gä enr) -1 
La 
À, 


Le În 
i emal(r) -1 
Li 
À À 
FAT DES 


T n=1,00 În 
i enf () -1 
rl 
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Simplifions ces expressions si cela est possible : 


(k | Ê ] | 
=| RENE) -l1 
L La 7 2 à, 
Tar 0) = F | (2 Sin(,0)= 27 Y A0 ( z) -| z) 


EEA În T o (2n+1 
7 (n+ l(i) e | CE, 
rl 


GLORI m 
roay a (2) Sn = F lÉ 2) {2)) 








r 








To (2n+1 


À, i 2, 
Eo 
4 > Sin(2,0) _ 27, 5 La r La r 
T n=0,% 2n+1 T n=0,% l À, 
(2n+1 (la) -1 
La 
2 qn Sin(4,0) ra DAE ON (2 Ÿ TE A 
in r ın r 
T ,0)=T. z 2— =| r E e n pha 
a2(7,0) 7 2 2n+1 | (=) (=) de E 2; 2n+1 (2) (a) | 


À À 
21, Sin(2,0)|{ r Pac 
T., (r,0) = T, — = 2 +| 
a= T À 2n+1 (2) (a) | 


Avec cette simplification vérifions au moins que la condition de continuité est encore valide, 
comme auto-contrôle de la validité des calculs : 


À 
DT en SALON 2 CE N” 

LOTO ME) en n7)} j| ta 
o C0) T 2 naa] (=) 
À 
mT emna NT 

T f ,0 = 0 n l= rl 
0 pm (| 


2 
i T 2 Sin(A,0 RE 
a20) = Tol) ET 7 E h- à | 








T2(r,0)=T, 




















nDo 2n+1 
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Admettons que 8,=n, alors il vient : 
À,=2n+1 








Lors 
2T. L La r fr 
To (r, 0) = =— D, Sin((2n+1)9) pourreQ,={0<r<1,,0<0<0,} 


T n= ,00 n 
i emf () -1 
rl 


À z 
2 (2) fg (2 
21 le r l, r 
To2(r,0)=T + — Sin((2n+1)0) pour r eQ, =$, <r<1,,0<0<0,) 


T° y=0,0 În 
Í Qn+1)| t2] -1 
La 


E 
rl + r2 2 
2T. G | G | | I E 
Taosa l a) Sin((2n+1)9) pourreQ,={l,,<r<w,0<0<0,} 


T 1,0 a 
en[r) 1 
rl 


Sin((2n + 15) = Cos(nx)=(-1)" 
(= v£) + (la) EH 2! 2n+1 
T 2T, r2 rl r 

Talr, 3) z 2 I 2n+1 | | 

(2n+1) (la) -1 

La 

f ” 2n+1 2 Qn+1 a $ 2n+1 L, 2n+1 

ou (2) DÉC | 


T 2T. 
To, )=T, + ° > 2n+1 
2 T n=0,% l 
(2n+ fe) _ j 
La 


I 2n+1 I 2n+l 
Ls l (+ | 
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Avec les expressions simplifiées : 


E 0) r 
rro- man | t | | 


J 
mee 5 CE a((2 j: 
z 








W aT (2n+1 


se Sin(,0) (a) 2276 Sin(2,0) Les 
ape) Re 2 2n+1 2- l J-( f hz T 2 2n+1 E LÉ | 


n=0,00 








En posant =n, on obtient : 


si (( (2 A 1) 2n+1 2n+1l in (À, of | | 
in((2n Sin 
0) = T 0)= 
AE > (2n +1) G | [2 -) | Er (2n +1) (2) Où 
2n+1 Qn+1 2n+1 2n+1 
2 = Sin((2n+1)0) r É | 2T, Sin(4,0)|{ r É | 
T(r,0)=T 2-|—| -|[# sh- : He 
o(,0)=T, = > De D 7 Re z DE b y 
Pour 0= - > Sin((2n + 1)8) =(-1)" 
T T 

Iar) E > 5 (D 5 2n+1 (r 2n+1 Tala) E 2 3 Cp” fl i (af 

T, Re (2n+ 1) l LS T n A (2n+ 1) r r 

2n+1 2n+1 

TaD 2y DS r (2) 

T, T aoo 2h ll ES r 


—] ñ 
Or + : ) l x"! est le développement de Taylor dela fonction Arctan(x) 
n=0,0 n+ 


























T 
Tass 


) 
E 2f ara =) z arc) pour re Q, = {0 <r<la0<0< 6,} 
TT rl r2 





To 





Ta2(r,0) _,_2 retan = + Aretan | pour reQ, ={, <r<la0<0<0, } 
Ty m [> r 


T 
Tr 


) 
—2 = 2 are 2) = areunf =) pour r EQ, = 1% <r<o,0<0< 6,} 
T T r r 


Les trois formes de la fonction sont identiques puisque la fonction Arctan, respecte la formule : 


1 
VxeR, Arctan(x)+ Arcta +) = > 
x 
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La fonction possède sur l'axe y, un maximum définie par : 


arf arc) 0 | l = l l > 
A S EN HE 
2 2 








ST 
~ 
ʻ 
N 
~ 
ah 
N 
+ 
x 
~ 
6 
N 
© 
+ 
oe 
~ 
N 
REET AN 
| = 
| 
| — 
SE 
Il 
~ 
x 
N 
| 
S— 
sy 
N 
| 
~ 
a 
z 
N 
`~ 
| 
LE 
© 


Problèmes aux limites de Laplace - secteur de disque en coordonnées polaires — p 153 


Exemple: Secteur de disque plein (r, 0) soumis à des conditions aux limites mixtes de Robin 
homogènes en Ÿ 


Pour le problème général : 
1 0 ÔT(r,0) 1 O’T(r,0) 

r + E 
r Or or ro 3o 
T(r,0)., fini 
T(r,0),.., =) 


asTo(r,0)+ B:T(r,0)),., =0 


0 





aTa (r,0)+ BTO), =0 


La solution est à rechercher sous la forme d'une série de forme générale : 
T(r,0)=e+ X r” (a,Cos(à,(0 -0,))+b,Sin(à,(6 -0,)) 


n=1,00 
Ce qui revient à appréhender des conditions aux limites polaires de la forme : 
aT, (r,0)+ B:T r0) =0 
a,T;(r,0)+ BTO) aa =0 
Notons au passage que ces deux conditions aux limites impliquent que e=0 dans la solution 
recherchée. D'après des calculs déjà réalisés dans le cadre de problèmes aux limites axi- 
symétriques, homogènes sur l'axe z. 
Fonctions propres polaires = ® (0) = a,4,Cos(4,(0 -0,))- B,Sin(2, (8 -0,)) 


T(r,0)= $ a,r*®,(0) 
n=1,00 $ 

Pour assurer la deuxième condition aux limites polaire homogène de Robin, il vient : 
À, valeurs propres solution de Sin(4,(0, -0, Ma aa 4 tB B,)= ,Cos(1, (6, -8, DIAA — p,a 1) 
Dans le cas des fonctions sinusoïdales, les normes se calculent ainsi : 
eo = 

1 (6, = X(,2,Cos(2, (6, —6, )- B,Sin(2, (6, —6, )) + (a, Sin(2, (0, -0, ))+ B,Cos(1, (6, —6, D) 
= 1 | 

2|+ z (&:4,Cos(à, (6, —6, )- B,Sin(à, (6, —6, )Xa,2,Sin(2, (6, —6, ))+ B,Cos(1, (6, —6, ))- a8; 


n 


La condition aux limites en r=l, donne : 
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[40 f,(@)&, (6) 


ST OE 


= = QÀ C -p;S, 
n=l,% Ti ®, O S CRO oi 


n n 


CRO 














n n 


0, 
B, = fao Jo (0) D,(8)=a;2 C —BS, 
0 


0, -0 


c=f dO f,(@)Cos(,(6-8,) ou C,= fdo f,(0+6,) Cos(2,6) 


0z “6 


S, = í dO f,(0) Sin(2,(0-8)) ou S,= fdo f,(0+0,) Sin(2,0) 
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Ce qui donne la solution suivante : 
2 
|e, O| = 


16, -0 N(æ:4,Cos(2, (0, -0,))- B:Sin(%,(0, -0,)} +(a,2,5in(2, (8, -0,))+ B:Cos(à, (0, -0,)F) 
E 2|+ L (aA, Cos(A, (0, — 6, )) T B:Sin(2, (0, -0, DXa;2,Sin(2, (0, -0, )) F B:Cos(å, (0, -0, )) —a;p; 





0, -0 


C, = Í d8 f (0) Cos(4,(0-0,) ou C,= Í dO f, (0 +0,) Cos(2,0) 


0> -0 


= [40 f0) Sin, (0-0) ou S,= [d0 f,(0+0,) Sin(2,0) 





rooy- J BASSE) ascosa (0-0))-Bsa, (0-8) 


Fr 
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Exemple : Secteur de disque plein (r 9) soumis à des conditions aux limites mixtes de Robin 
homogènes en Ÿ 


Pour le problème général : 
1 8 ÔT(r,0) 1 Ô°T(r,0) 

r + = 
r ôr ôr ro 80 
T(r,0) „fini 
aT; (r,0)+ BT r, _, = f,(0) 


aT, (r,0)+ BTO) =0 


0 








aT, (r,0)+ BTC, O) p =0 


D'après des calculs déjà réalisés dans le cadre du problèmes aux limites précédent, les deux 
dernières conditions aux limites polaires donnent : 

Fonctions propres polaires = ® (0) = a,4,Cos(4,(0 -0,))- B;Sin(à,(0 -0,)) 

À, valeurs propres solution de Sin(A, (0, —0, Maaa a +b: B,)= ,Cos(1, (0, — 06, Ba; — Ba) 
T(r,0)= Y'a,r"®,(8) 


n=1,00 
Dans le cas des fonctions sinusoïdales, les normes se calculent ainsi : 
2 
(e, O = 


T0, -aX(a:2,cos(2,(8, -8,)-B,5in(2, (0, -0,)) + (a32,Sin(à, (0, -0,))+ B.Cos(z, (0, -0,)}) 
f 2|+ ; (a;À,Cos(1, (6, -0, )- B,Sin(2, (6, 4 JXa2,Sin(2, (6, —6, ) + B,Cos(à, (6, —6, ))- ap; 


n 


La condition aux limites en r=l, donne la solution suivante :: 
aT, 0) + pT (r,0)| _ sal" "+8 Yael" + p, 9, O= AO) 




















n=1,00 

0, 
fdo f0) D, 0) 

= Q = B, = À, C, = b; S, 

"Rat + pa NOO (34,1, + BL EO (4, + Ba, EO 
0, 0, CE 
B, = | d0 f,(@)®,(6)=a2,C,-BS, C,= [40 f,(0)Cos(,(0-8) ou C,= [d0 f,(0+0,) Cos(4,0) 
0, CA 0 


S, = Î dO f(O) Sin(2,(0-8) ou S,= [do f,(0+0,) Sin(2,0) 


0, 





l 
T 2| +> (a4, Cos(4, (0, -0,))-B;Sinl2, (0; -8))Ke,2,5in(2, (0, -8))+ 2:Cos(4,(0 -0,))- a8; 
(a,2,C,=B,S,) 
T ,0 = n n n 3 
na Lars B.1" EO] 





r#(a,À,Cos(2,(0 -6,))- B,Sin(2,(8 -0,))) 
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Exemple : Secteur de disque annulaire (r,0) soumis à des conditions aux limites de Dirichlet 
homogènes en 0 


Occupons maintenant d'un secteur de disque annulaire soumis à des conditions de Dirichlet sur 
chacune de ces frontières. Sans restreindre le problème on peut fixer le rayon intérieur à 1 et le 
rayon extérieur à k. On place également le secteur annulaire le long de l'axe des abscisses. 


2 
1 ô TC) 1 TOA) 


0 
r Or Or r” 00? 








T(r,0)|_. =0 
T(r,0)|., =1 
T(r,0),,, =0 
T(r,0),. =0 


On rappelle que les solutions de l'équation séparée en fr,0) sont de la forme 


T(r,0)=e+ fin(r)+ X (c,r* +d,r * \a,Cos(2,0)+b,Sin(2,0)) 


n=1,00 


La condition de nullité en 0=0 donne an=0, il n'y a plus de condition de finitude de la solution et la 


condition en r=1 donne ® Faig De la configuration des conditions aux limites en r=1 et 0=0, 


réunion de deux conditions nulles sur les tranches du secteur, on déduit facilement que la valeur 
de e est nulle . La solution prend la forme plus simple : 


T(,0)= fin(r)+ Xb, (> -r> ]Sin(2,0) 
n=1,00 
a 


Appliquons maintenant la conditions aux limites en =a, 


Il suffit finalement d'appliquer la condition aux limites en r=k, pour trouver la solution du 
problème, en utilisant les relations d'orthogonalité des fonctions propres ©(0,A,)» SUr l'intervalle 


0e [6,,8,]- [0.27] soit: 


| do 6(.2,8(0.2,,)- 0 si n + m 
O(0,2,)=Sin(2,0) ou @(60,À,)=Cos(2,0) 
1-Cos(20)_ 1 f s 





Sina O| = p dO Sin? (2,0) = p d0 





2 An 2 4 1,0 
> Amba 
|Cos(2,0)[ = j dO Cos? (1,0) = À go + Cos(28) _ 1 [8 _ Sin(26) 
0 0 2 À, 2 4 1,6 


Problèmes aux limites de Laplace - secteur de disque en coordonnées polaires — p 158 


D D, (k> -k> JSin(2,a) =1 


n=1,00 








S 7 1na ; nī 
[ d0 Sin? (4,0) = 1 É seo _ É seo _Aa_a 
0 À, | 2 4 j Al2 4 73i -3 


f . 1-Cos(i,a) 
2f 40 Sin(2,0) fa | cent 


230 D n = 
fe alk k>) alk- ™)] Aak -k™) 
4 





>b,,=0 





(2n + 1)}r 4 (r> — pr }sin(à 0) 
À => , T. ,0 _ n 
n à 28,1 (r,0) a À (n+1fk” k”) 


Pour le problème : 

13 ÔT(r,0) 1 &T(r,8) 
r + = 

r ôr or r ô’ 

T(r,0) _ =1 

T(r,0),., =0 

T(r,0),.. =0 

T(r,0)|,_, =0 


On trouverait de même : 


T(r,0)= 2 IA - e J smao) 


CLenr=k> Ý -b,(k* -k> )Sin(4,a) =1 


n=1,% 





0 











fa Sin? (4,0) = 5 
2f d0 Sin(2,0 
f Í in(4,0) i 21-9") opasi 4 


== 0 n Z 
n` alk” -k ) un 2, al" xx) 2n+l — 2, ak -k ) 
(2n+1)r 


a | nom-s g CEM pan 


no (2n+ ie fn ) 





Problèmes aux limites de Laplace - secteur de disque en coordonnées polaires — p 159 


Ce qui donne pour le problème, par principe de superposition : 


2 
1 0 .6T(,6), 1 0!T(,6) 


= (0 
r Or Or r” 00? 





T(r,0),_. =7 
T(r,0),., =T 
T(r,0),.,=0 
T(r,0),_, =0 





T(r,0)= DE (72) - k i }sin,0) 
L3(7,0) i D1(,0)+T Ts (,0) 


À —4, i "n k/ |" ; 
EEEF 4 D (r "n y " JSin(2,0) T 4 5 (Al VA] }sin,0) 


T 4. Onte” + ra. (n+1K-K*) 


nets a aCi r28) 


r (2n +1k* 74 








Là encore on ne peut trouver de solution par la méthode de séparation des variables lorsque les 
conditions aux limites cesse d'être homogène en 06. 


En conclusion de cet exemple on peut aborder le problème plus général : 


2 
1 0 .OT(,6), 1 0?T(r,6) 





r Or Or r” 00? 
T(r,0) Pha 

TE, 0), = f(0) 

T(,0), = 0 

T(r,0) Fe 





0; “6, 


b = Í d0 f,(0) Sin(4,(0-8) ou b,= Í dO f,(0 +0.) Sin(4,0) 
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Exemple : Secteur de disque plein (r,0) soumis à des conditions aux limites mixtes de Dirichlet, 
Neumann et Robin homogènes en 6 


Occupons nous maintenant d'un secteur de disque plein soumis à des conditions de Dirichlet sur 
chacune de ces frontières. Sans restreindre le problème on peut fixer le rayon du disque à 1. 


13 ÔT(r,0) 1 Ô°T(r,0) 

r + ` 
r Or Or r” 00? 
T(r,0)., fini 
T(r,0) = 1 
T(r,0),., = (0 


Ti(r,0) _ =0 


0 











La solution est de la forme : 


T(r,0)=e+ fin(r)+ X (c,r* +d,r * Ja, Cos(2,0)+b,Sin(2,0)) 


n=1,00 


La condition de nullité en 6=0 donne an=0, la condition de finitude de la solution donne dn=0 et 
f=0. De la configuration de finitude en r=0 et de la condition nulle sur une des tranches on en 
déduit que la valeur de la solution en r=0 est nulle également, donc e=0. La solution prend la 
forme plus simple : 


T(r,0)= $ b,r” Sin(2,0) 


n=1,00 
Appliquons maintenant la conditions aux limites en O=a 


(2n+1}r 
2a 


La condition en r=1 donne : 
$ b,Sin(4,0)=1 


n=1,00 


Cos(A,a)=0= À, = 


(2n+#1) 
1 f eo |" zA f 5 1a 
À. |2 4 n A, L2 4 0 24 


n n n 








[ 46 Sin°(2,0)= = 
J 2 
2[ 46 Sin(,0) f emam) 
b je md 0 ps À, pr 2 
g a a A,a 
(2n+1)r 4 = r”Sin(4,0) 
= , I,,(r,0)= r 
5 2a 2408) T > (2n+1) 








Problèmes aux limites de Laplace - secteur de disque en coordonnées polaires — p 161 


Et pour le problème 

2 
1 ə ,2T00) | 1 ð T(r,0) o 
r Or or r ô’ 
T(r,0) „fini 
T(r,0)|. =1 








T(r,0),_ =0 

Sur la solution de la forme (on exprime de la même manière les fonctions propres avec le 
changement de variable 0->a-0) 

T(r,0)=e+ fin(r)+ X (c,r* +d,r * Ja, Cos(à, (a -6))+b,Sin(à, (a -0)) 


n=1,00 


La condition de finitude de la solution donne dn=0 et f=0. De la configuration de finitude en r=0 et 
de la condition nulle sur une des tranches on en déduit que la valeur de la solution en r=0 est nulle 
également, donc e=0. La condition de nullité en 0=a donne an=0. La solution prend la forme plus 
simple : 
T(r,0)= X b, r” Sin(A, (a -0) 

n=l . Au passage vérifions les conditions d'orthogonalité des fonctions 
sur [0,a] qui demeure équivalente en utilisant le même changement de variable 0->a-0 : 
g'=a-ÿ—= 


[ d0 Sin(À, (a — 9))Sin(À,, (a - 9)) = — [as Sin(À,9')Sin(À, 9") 


= [ d9' Sin(À,9')Sin(À, 9") = 0 si n# m 





$ Ama 
]Sin(4, (a 9) = [ d0 Sin” (À, (a -9)) = -fas Sin? (A9) — - É seo) 


2 4 


0 
Appliquons maintenant la condition aux limites en 0=0 


2n+1 
Cos(l,a)=0= À, = (@n+1)r 
2a 
et en utilisant ces relations d'orthogonalité, la condition en r=1 donne : 
D b,Sin(2, (a -0)) =1 


n=1,00 
[ d0 Sin? (2, (a -0)) = | d0 Sin? (2,0) = 5 
0 0 


2f d0 Sin(,(a -0) 2f d0 Sin(,0) 





b, = ————— = < -—— 
a a 2,a 
De (2n+1)r | Tt.0)- 4 5 r'"Sin(,(a -0) 
2a Mn (2n+1) 
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Le problème : 


2: 
1 0 ,0T(r,0) | s 0 BA 
r Or or r 00 
T(r,0)|., fini 











T(r,0),. =1 
Tr, =0 
T, rO =0 


a une solution triviale T,„,(r,0)=1 


Toutes les solutions aux problèmes généraux avec une condition aux limites fa(9) sont les suivantes: 


Pour le problème : 
10 OT(r,0) 1 Ô°T(r,0) 

r + = 
r Cr ôr r 00 
T(r,0) fini 
T(r,0),.., = A0) 
TO), =0 


T(r,0) =0 


0=0, 





0 








0 -0 


b, = Í dO f,(6) Cos(4,(0-0,) ou b,= [d0 f,(0+6,) Cos(2,0) 





PE TED =T wE) b,Cos(4,(0-0,) 


r 


Ou encore le problème suivant : 
1 0 OT(r,0) 1 Ô°T(r,0) 

r + = 
r Cr ôr ro 00 
T(r,0)| fini 
T(r,0),., = f0) 
T(r,0),., = 0 





0 





T(r.0),., =0 


® 6 


b, = Í dO f(O) Cos(à,(0-06,) ou b,= [d6 f,(0,-0) Cos(4,0) 


_Qn+Dz ME rŸ Í 
A= a] f reo- lE] b,Cos(4,(0-0,) 


r 


Ou encore le problème suivant (ajouter la fonction de valeur propre nulle !): 
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1 0 OT(r,0) 1 Ô°T(r,0) 
r + = 

r Cr Or r” 00? 

T(r,0)., fini 

T(r,0),., =A) 


T, (r,0) 


0 











= 0 
0=6, 


T(r.0),., =0 
6,-6, 


f= Í d0 f,(0) b,= í dO f,(6) Cos(à,(0-0,) ou b,= fdo fa(0+0,) Cos(2,0) 








O AT o f 2 E B 
À, = l LTEN zoz b,Cos(4,(0-0,) 


r 
Passons au problème suivant avec des conditions aux limites mixtes de Robin 


2 
1 ô TC) 1 TCA) 


0 
r Or or r 00? 











T(r,0) „fini 
T(r,0). =1 
T(r,0),., = 0 


T,(r,0)+AT(r,0)). =0 
sur la solution recherchée : 
T(r,0)=e+ finr)+ X hc,r* +d,r™ (a, Cos(4,0)+b,Sin(4,0)) 


n=1,00 


La condition de nullité en 6=0 donne an=0, la condition de finitude de la solution donne dn=0 et 
f=0. De la configuration de finitude en r=0 et de la condition nulle sur une des tranches on en 
déduit que la valeur de la solution en r=0 est nulle également, donc e=0. La solution prend la 


forme plus simple : 
T(r,0)= $ b,r*™ Sin(2,0) 


n=1,00 


Appliquons maintenant la conditions aux limites en =a, cela donne 


À,Cos(2,a)+h Sin(,a) = 0 


C'est l'équation qui définit les valeurs propres du problème de Sturm-Liouville. Par exemple en 


prenant a=n/4 et h=20, voici les trente premières valeurs propres : 
3.7632, 7.54092, 11.3431, 15.1736, 19.0316, 22.9137, 26.8159, 
30.7345, 34.6663, 38.6086, 42.5593, 46.517, 50.4803, 54.4482, 58.42, 
62.3949, 66.3726, 70.3527, 74.3346, 78.3183, 82.3035, 86.29, 90.2776, 
94.2662, 98.2557, 102.246, 106.237, 110.229, 114.221, 118.213 
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La condition en r=1 donne : 
$ b,Sin(4,0)=1 


n=1,00 


fao si (2,0) = 12- Le =- [2e eaa). fa sasaa 








2 4 2 4 2 


À, n 


nanam) 
ie à, ___ 2(1-Cos(4,a)) 


a DOLA) (q MCAD) 3 a Sin22,a) 
24, 2A, 24, 


À, solutionde À he a)+h Sin(à,a)=0, 


s 25 "(1- Cos(i,a))Sin(2,0) 
Tas E | D 
AE ei a 
2A, 


2 Í d0 Sin(,0) | 





Et pour le problème 

1 0 ÔT(r,0) 1 OT(r,0) 
r + 

r Cr or r 00 

T(r,0)., fini 

T(r,0) =1 

T,(r,0)-AT(r,0) = 

T(r,0), =0 


O=a 


=0 








La solution est la suivante en effectuant le changement de variable 0->a-8 puisque pour 
T(r,0)= Y°b,r”"Sin(2,(a -0) 
ee , la condition aux limites en 0=0, donne une équation aux valeurs 


propres identique — À,Cos(2,a)-h Sin(2,a)=0< 2,Cos(2,a)+h Sin(2,a)=0 


À, solution de À ui a)+h Sin(À,a)=0, 


T, 0,8) = ST "(1— Cos(À,aæ))Sin(2, (a —-0)) 


n=1,% Ilas Sin(22,a) 
22, 
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Pour le problème plus général : 
2 
Lo OT(r,0) pl © 1050) -0 
r Cr ôr Ñ 00 
T(r,0)| fini 
T(r,0),.., =) 
Tr, O) =0 
T,(r,0)+AT(r,0) =0 


La solution se développe en série de fonctions propres : 








2 fao Jfa (0) Sin(2,0) a 
b. = — b = | d0 0) Sin(À 0 
n Eo > D, | Jo (0) Sin(å,0) 
2A, 


À, solution de  À,Cos(i,a)+h Sin(2,a) = 0, 


T 4(r,0) =2 > | b, Sin(2,0) f n 





n=l% Sin(22,a) 
Q — 
24, 


et le problème complémentaire : 
1 0 OT(r,0) 1 Ô°T(r,0) 

À t- ere 
r Cr ôr à 00 
T(r,0) fini 
T, 0), =A) 
T,(r,0)-AT(r,0) =0 
T(r,0), =0 


O=a 





0 





La solution se développe en série de fonctions propres : 


z fao f,(0) Sin(A, la -0) à 


b = b = |d0 0) Sin(À (a —-0 
n E >D, | Jol ) Sin( „(a ) 
2À, 


À, solution de  À,Cos(i,a)+h Sin(À,a)=0, 


b, Sin(,(a-6) | r j 
L 








Tosa(r:0) E 2 TE Sin(22,@) 
2, 
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Pour le problème mixte suivant : 


1 0 ÔT(r,0) 1 O’T(r,0) 
r + = 

r ôr or r? 80’ 

T(r,0),., fini 

T(r,0)., =1 


Ti(r.0),., =0 


0 














T(r,0)+ AT(r,0) =0 


sur la solution recherchée : 
T(r,0)=e+ fin(r)+ X (e,r +d,r™ \a,Cos(2,0)+b,Sin(2,0)) 

n=1,00 
La condition de nullité de la dérivé première en 0=0 donne bn=0, la condition de finitude de la 
solution donne dn=0 et f=0. De la configuration de finitude en r=0 et de la condition nulle sur une 
des tranches on peut dans un premier temps faire l'hypothèse que la valeur de la solution en r=0 
est nulle également, soit e=0. La solution prend la forme plus simple : 
T(r,0)= $ b,r™ Cos(4,0) 


n=1,00 


Appliquons maintenant la conditions aux limites en O=a, cela donne 


—4,Sin(4,a) +h Cos(,a)=0 C'est l'équation qui définit les valeurs propres du problème de 


Sturm-Liouville. Par exemple en prenant a=n/4 et h=20, voici les trente premières valeurs 
propres : 

1.88063, 5.64948, 9.43858, 13.2548, 17.0993, 20.9699, 24.8625, 

28.7734, 32.6989, 36.6362, 40.583, 44.5374, 48.498, 52.4637, 56.4336, 

60.4071, 64.3835, 68.3624, 72.3434, 76.3263, 80.3108, 84.2966, 

88.2837, 92.2718, 96.2608, 100.251, 104.241, 108.233, 112.225, 116.217 


La condition en r=1 donne : 
$ b,Cos(4,0)=1 


n=1,% 


œ š Aa : | 
f 46 Cos”(1,8) = 1 É . Se) sd [ze E furent) _1f,,, Sin(22,a) 
0 À, L2 4 0 A,\ 2 4 2 21 


n n n 








2f d0 Cos(1,0) daaa) 


n 


E E 2Sin(å a) 
PE Sin(22,a) Pe Sin(22,a) Ahaa Sin(22,a) 
22, 2A, 2A, 
À, solution de A Sin(2,a) = h Cos(2,a), 


r” Sin(å,æ)Cos (4,0 
Daey, ( n ) ( n ) 


n=1,00 À, HSE Sin(2,c) 
22 


n 
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Et pour le problème mixte 


10 OT(r,0) 1 Ô°T(r,0) 

r + = 
r Or or r? 00? 
T(r,0)| „fini 
T(r,0) 


T,(r,0)-AT(r,0) = 0 


0 








=] 
PSI 


T, (r,0) =0 


en effectuant le changement de variable ĝ->a-Ô dans la fonction précédente on obtient la solution 
recherchée : 


À, solution de À,Sin(à,@)=h Cos(A,@), 
À, © A 
T(r,0)=2Y Sin(A,æ)Cos(2, (a -0)) 


n=1,00 À a+ Sin(22,@) 
k 2A 


n 





Pour le problème plus général avec ces conditions aux limites mixtes Neumann et Robin : 





1 0 OT(r,0) 1 Ô°T(r,0) 
r + = 

r Cr Or ro 00? 

T(r,0)| fini 

T(r,0),, = f.(6) 


= 0 
0=0 


0 





T,(r,0) 





T(r,0)+ hT(r,0) =0 


La solution se développe en série de fonctions propres 
À, solution de A Sin(2,a) = h Cos(A,a), 


CL. Yb,l"Cos(2,0)= f,(0) 


n=1,0 








fao Cos (0) =2 a + Fees) 
J 2 2, 
2[ 40 f, (6) Cos(,0) N 
p- E, 
: | sx) =: j Ja (0) Cos(2,0) 
RS, à 





À 

b, Cos(4,0) (rY 

Te D Sin(22,@) É 
Par [ST 


2% 
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Et pour le problème complémentaire avec ces conditions aux limites mixtes Neumann et Robin: 


16 ÔT(r,0) 1 @T(r,0) 
r + = 

r Or or r” 00° 

T(r,0)| fini 

T(r,0),., = f.(0) 

T, (r,0) -hT (r0) =0 


0 





T;(r,0)|_,=0 
Nous obtenons la série de fonctions propres 
À, solution de À Sin(A,&)= h Cos(à,æ), 


b, 5 fao fa) Cos(1, (a 0) 
T5,6(7,0) = 2 > b, Cos(A, (a z) n 


n=1,00 a+ Sin(22,@) 
2A, 





Exemple : Secteur de disque annulaire (r,9) soumis à des conditions aux limites mixtes de 
Dirichlet, Neumann et Robin homogènes en 0 


Occupons nous maintenant de problème similaire d'un secteur de disque annulaire soumis à des 
conditions de Dirichlet sur chacune de ces frontières. Sans restreindre le problème on peut fixer le 
rayon intérieur à 1 et le rayon extérieur à k. 


2 
10 ,2T(r,0) | 10 T(r,0) o 
r Or Or r” 80’ 
TO) =0 


T(r,0)., =1 
T(r,0),., =0 
T(r,0). =0 





La solution est de la forme : 


T(r,0)=e+ finr)+ X (c,r* +d,r™ (a, Cos(4,0)+ b,Sin(1,0)) 


n=1,00 


La condition de nullité en 0=0 donne an=0, f=0, la condition en r=1 donne cn=-dn pour n>0 et e=0 
pour n=0. La solution prend la forme plus simple : 


T(r,0)= $ b, (r> =r% ]Sin(4,0) 


n=1,00 
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Appliquons maintenant la conditions aux limites en O=a 
(2n+1}r 
2a 
La condition en r=k donne : 
D D, (k> — 4% ]Sin(2,0) =1 


n=1,00 


Cos(1,a)=0= À, = 


fao Sin?(2,0)=© 
J 2 
2f d0 Sin(2,0) 





b,=— = 
7 alke-f#) Aak--##) 


_(en+iÿr L4 y (è -r> Sma) 
a 2a + Rte m 2 (2n +1\(k* xx) | 





Et pour le problème 





2 
10 ,2T(r,0) | 10 T(r,0) o 
r Or Or ro 80’ 





T(r,0)|__=0 
T(r,0)., =1 

T, CRD =0 
T(r,0),. =0 


Sur la solution de la forme (on exprime de la même manière les fonctions propres avec le 
changement de variable 0->a-0) 
T(r,0)=e+ f In(r)+ Sr” +d r^” \a,Cos(2, («-0))+b,Sin(, (a -0)) 


n=1,0 


La condition de nullité en =a donne an=0, f=0, la condition en r=1 donne cn=-dn pour n>0 et e=0 
pour n=0. La solution prend la forme plus simple : 


T(r,0)= Y b,(r* =r% ]Sin(à, (a -0) 


n=1,00 
Appliquons maintenant la condition aux limites en 0=0 
(2n+1}r 


Cos(i,a)=0—= À, = 
2a 


Problèmes aux limites de Laplace - secteur de disque en coordonnées polaires — p 170 


et en utilisant ces relations d'orthogonalité, la condition en r=k donne : 
D b, (k> -K> Sina, (a -0) =1 


n=1,00 
[48 Sin? (2, (a -0) = [49 Sin? (2,0) => 
0 0 


2f d0 Sin(A,(a -0) 2f dO Sin(,0) 











b = — = - 
Earl dé) a Tr 
_(2n+1}r C4 g (ré -r> Sina, la -0) 
a 2a > To208)= m 2 (2n+ 14> -k^ ) 
Le problème : 


1 0 ÔT(r,0) 1 @T(r,0) 

r + = 
r ôr or ro 80’ 
T(r,0)., =0 


T(r,0) 
T, (r,0) 


0 








r=k =] 


= 0 
0=0 





T ro) =0 


a une solution triviale T;,,(r,0)=1 


Toutes les solutions aux problèmes généraux avec une condition aux limites f 4{ 0) sont les 
suivantes : 


Pour le problème : 
1 ə ,2T0:0) | 1 &’T(r,0) B 


0 
r Or Or r” 00? 








T(r,0),., =0 
TEO), = fo (0) 
T(r.0),., =0 
T(r,8) pe, =0 





6 -0 


b, = Í dO f (0) Cos(4,(0-0,) ou b,= Í dO f (0 +0,) Cos(4,0) 
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Ou encore le problème suivant : 
13 ÔT(r,0) 1 Ô°T(r,0) 

r + 5 7 = 
r Or ôr r 00 











T(r,0)., =0 
TC), = h0) 
TO) p =0 


Tr), =0 


0, -0 


b, = Í d0 f (0) Cos(A,(8-0,) ou b,= Í d0 f,(8, -0) Cos(2,0) 


_ (n+l) 2 


aao) * ea 2 b,Cos(3,(0-0,) 


(0, -0,),5. TANE 
La L 
Ou encore le problème suivant (ajouter la fonction de valeur propre nulle !): 


13 ÔT(r,0) 1 &T(r,8) 
r + = 
r Or Or ro 80’ 








0 














T(r,6),., =0 
T(r,0),... = fo) 
TEO) = 
T, CO) p =0 


0, “6 


f = Í d0 f,(8) b,= Í dO f,(0)Cos(4,(0-0,) ou b,= Í d0 f (0 +0,) Cos(4,0) 











AR) 

nm f r 2 La r 

À = 5 T. ,0 = l b,C À, 0 0 

EO 2e) 1 ea AAA) Ee 
BI 


Si les conditions aux limites sont inversées entre l, et l, il suffit d'effectuer les changements : 
l:->1, et l,,->l, dans les expressions précédentes. 
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Passons au problème suivant avec des conditions aux limites mixtes de Robin 
2 

1 © 2700) | à Ô UE) Li 

for or r 00 

TO) =0 

TO =] 

T(r,0),., = 0 

Ti(r,0)+ AT(r,0),. = 0 

sur la solution recherchée : 

T(r,6)=e+ finr)+ X (c,r* +d,r™ (a, Cos(,0)+b,Sin(1,0)) 


n=1,00 





La condition de nullité en 6=0 donne an=0 et f=0, la condition en r=1 donne cn=-dn pour n>0 et 
e=0 pour n=0. La solution prend la forme plus simple : 
T(r,0)= D b,(r* =r% ]Sin(2,0) 


n=1,00 


Appliquons maintenant la conditions aux limites en @=a, cela donne À, Cos(,a)+ h Sin(,a) = Le 


C'est l'équation qui définit les valeurs propres du problème de Sturm-Liouville. Par exemple en 
prenant a=n/4 et h=20, voici les trente premières valeurs propres : 

3.7632, 7.54092, 11.3431, 15.1736, 19.0316, 22.9137, 26.8159, 

30.7345, 34.6663, 38.6086, 42.5593, 46.517, 50.4803, 54.4482, 58.42, 

62.3949, 66.3726, 70.3527, 74.3346, 78.3183, 82.3035, 86.29, 90.2776, 

94.2662, 98.2557, 102.246, 106.237, 110.229, 114.221, 118.213 


La condition en r=k donne : 
Yp, (k> -k> )Sin(4,0) =1 


n=1,00 


a x Ana à Y 
fao Sin? (4,0) = l É seo) z l (as sea). 1 m Sin(22,@) 
A À, [2 4 ò À 2 4 2 2 


ae 








n 


À 


n 


] 2 
2 Í dO Sin(2,0) | | eo 


b = 0 = = 
, [a OLA es p) (a OAA es p>) ala OAA fee p>) 
24, 21, 2À, 


À, solution de  À,Cos(2,a)+h Sin(A,a) = 0, 

r —r™ 1 Cos(A.a))Sin(2,0 
2609-25 | 1 T ` )) ( n ) 
n=1,00 ala a ın „æ je K>) 




















2A, 
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Et pour le problème 

2 
1 © 2700), 1 Ô T(r,0) o 
r or or r ô’ 
T(r,0),._ =0 


T(r,0)., =1 
T(r,0) —AT(r,0). =0 
T(r,0),. =0 








0= 


La solution est la suivante en effectuant le changement de variable 0->a-8 puisque pour 


T(r,0)= $ b, (r* -r™ Sina, (a -0) 
nake , la condition aux limites en 0=0, donne une équation aux 


valeurs propres identique À, Cos(,a)-h Sin(å,a)=0 


À, solution de  À,Cos(2,a)+h Sin(2,a) = 0, 

r® —r™ (1- Cos(A,a))Sin(A, (a -0 
rat =2F | k- a 2 G,(a-0) 
n=1,0 À, (o £ POLA g =- k>) 





24 


n 


Et pour le problème mixte suivant : 





2 
1 ô TC) 1 TCA) 


0 
r Or or ro 80’ 








T(r,0).. =0 
T(r,0)., =1 
T,(r,0),.=0 





T,(r,0)+AT(r,0). =0 
sur la solution recherchée : 
T(r,0)=e+ fin(r)+ X (e,r +d,r™ \a,Cos(2,0)+b,Sin(2,0)) 
n=1,00 
La condition de nullité de la dérivé première en 0=0 donne bn=0, la condition en r=1 donne cn=- 
dn pour n>0 et e=0 pour n=0. La solution prend la forme plus simple : 


T(r,0)= fin(r)+ Xb, (r* -r JCos(2,0) 


n=1,00 


Appliquons maintenant la conditions aux limites en  @=a, cela donne f=0 pour n=0 et 


A Sin(,a)+h Cos(i,a) = 0 pour n>0. C'est l'équation qui définit les valeurs propres du 
problème de Sturm-Liouville. Par exemple en prenant a=n/4 et h=20, voici les trente premières 
valeurs propres : 

1.88063, 5.64948, 9.43858, 13.2548, 17.0993, 20.9699, 24.8625, 

28.7734, 32.6989, 36.6362, 40.583, 44.5374, 48.498, 52.4637, 56.4336, 
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60.4071, 64.3835, 68.3624, 72.3434, 76.3263, 80.3108, 84.2966, 
88.2837, 92.2718, 96.2608, 100.251, 104.241, 108.233, 112.225, 116.217 


La condition en r=1 donne : 
Yp, (k> -k> JCos(4,0)=1 


n=1,0 


a R Aa A R 
fao Cos’ (4,0) = l É r Sen P (2e į ee) A (a { Sga) 
0 




















AE AE 4 2 22, 
2[ d0 Cos(2,0) {2a AA 
b, = 0 = 'n = in „æ 
ad Sin(2,c) (k> k> ) gi Sin(22,a) (r> k ) À ği Sin(2,a) (x p> ) 
TR À. FY 


À, solution de h Cos(A,a)=2,Sin(2,a), 
r^ —r ” ]Sin(A,a)Cos(2,0 
To9(r,0)=2 0 a Cu 0) 


n=1o les Sin(22,@) (k> Lg) 
2 


n 





Et pour le problème mixte complémentaire 





2 
1 ô TC) 1 TCA) 


0 
r Or Or ro 80’ 








T,(r,0)-hT(r,0)| _ =0 
T,(r,0)_. =0 


en effectuant le changement de variable 8->a-8 dans la fonction précédente on obtient la solution 
recherchée : 


À, solution de h Cos(2,a) = À, Sin(A,a), 


An An |: E 
egar C a a (a -0) 
nié a [e + POA Vus LEA ) 


n 
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Les solutions aux problèmes généraux avec ces conditions aux limites mixtes Dirichlet, Neumann et 
Robin sont : 





2 
1 ô OT(r,0) 1 TE) 


0 
ror Or ro 00? 








T(r,0)  =0 
T(r,0),.,, = f9(0) 
Tr, Op =0 


Ti(r,0)+AT(r,0),., =0 


a =0,-—0, ; À, solution de  À,Cos(,a)+h Sin(A,a) = 0, 


bre Í d0 f (0) Sin(4,(0-0,) 


> BI F) 
; L, 
L5:10:0)=2 b, Sin(2,0) r r 


7 À, À 
21, La La 


et pour le problème complémentaire : 





2 
10 ôTC,0) l ð TED -p 
r Or Or r 00 
T(r,0)., =0 
TEO), = fo (0) 
T,(r,0)-AT(r,0) =0 








T(r, 0), = 0 


on a en effectuant le changement de variable Ü-Ü1->a+01-d soit 0-01->02-0 
a =0,-—0, ; À, solution de  À,Cos(,a)+h Sin(A,a) = 0, 


b, = | d0 f,(0) Sin(2,(. -0) 


Á i La * 
E 
Ton (r,80)=2 > b, Sin(2, (6, -0)) | r 


2T ala z gam) { A j Š I f) 
24, l L, 
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Pour le problème 
1 ð ,2T(r,0) | 1 O°T(r,0) S 


0 
ror Or r” 00? 








T(r,0).. =0 
T(r,0),., = f(0) 
TO) p =? 





T,(r,0)+ AT(r,0),., =0 
ond: 


a =0, -0 ; À, solution de h Cos(,a) = À,Sin(2,a), 


b, = [d0 f,(6) Cos(2,(8-0,) 


á r la i 
Moe 
Toi (r,0)=2 > b,Cos(2, (0 8,) | 7 k 


. À À, 
ar (agaa) à) { La 
2A, La La 
et pour le problème complémentaire : 
2 
1 ð , 2T(,0) , 1 OC T(r,0) i 
r ôr or ro 8o 
T(r,0) ; =0 
TEO =O 


To(r,0)—hT(r,0),., =0 














Tyr, O =0 


on a en effectuant le changement de variable 0-01->æ+01-0 soit 0-01->02-0 
a =0,—0, ; À, solution de h Cos(A,a)=2,Sin(2,a), 


b, = [d0 f,(6) Cos(à, (8. -6) 


4 à La i 
ao) 
L013(7,0) = 2 >. b,Cos(4, (6, 0) | 4 


= fes meaa) wy [ay 
22, La La 





Lorsque l'on inverse les conditions aux limites dur ln et l, il suffit de réaliser les changements : 
la -> l2et l2 -> ln pour obtenir les solutions des problèmes généraux. 
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Exemple : Secteur de disque annulaire (r,Ÿ) soumis à des conditions aux limites mixtes de Robin 
homogènes en Ü, inhomogènes en r 


Pour le problème général : 

1 0 ÔT(r,0) 1 O’T(r,0) 

— — r — + — ————— — 
r ôr or ro 3o 
T(,0),., = 0 

T(r,0),.. =A) 


asTo(r,0)+ B:T(r,0)),., =0 





oTo (r,0)+ BTO), =0 


D'après des calculs déjà réalisés sur les fonctions propres polaires dans la prise en compte des 
deux conditions aux limites : 
Fonctions propres polaires = ® (0) = a,4,Cos(4,(0 -0,))- B,Sin(2, (0 —0,)) 


À, valeurs propres solution de Sin(4,(0,—0, a, aa 4 +p; B,)= 2,Cos(1, (6, -0, )\B,a, - B.a,) 


r 3 La a 
roD- J ad) (2) ko 


Dans le cas des fonctions sinusoïdales, les normes se calculent ainsi : 
(D, @) = 

6-0, X(c,2,Cos(2, (8, -0,))- B,Sin(2, (0, -0,)} + (@,4,Sin(2, (0, -0,)+ B,Cos(, (8, -8,))) 
Tal — (a,4,Cos(4,(0, -0,))- B,Sin(2, (0, 0, ))\a,1,Sin(2, (8, -0,))+ B,Cos(2, (8, —-8,)))-a,p, 


n 


La condition aux limites en r=l,; donne : 











Pegy fado j,(6)®, (6) 
TOE ER 
FES GT 
— b, — a2,C, Bi, 
2 E Nh | mue | 
(a) er (ehe 


n n 


0, 
B, = | d0 f,(0) D, (0)=0,4,C, -B:S, 
0 


0> -0 


C, = Í dO f,(0) Cos(à,(0-0,)) ou C,= fdo f,(0+6,) Cos(2,6) 


0, -0 


S, = í dO f (0) Sin(4,(0-0,) ou S,= f dO f,(0 +0,) Sin(À,0) 


Problèmes aux limites de Laplace - secteur de disque en coordonnées polaires — p 178 


Ce qui donne la solution suivante : 


lo (0) = 





L, ((,4,Cos(4, (6, —6, ))- B,Sin(2, (6, —6, D) + (a,2,Sin(à, (6, —6, ) + B,Cos(4, (6, —6, )} ) 


1 
E 2|+ _ (asà, Cos(2, (0, -0, )) a B:Sin(à, (, -0, DXa;2,Sin(2, (0, -0, )) h B;Cos(A, (0, -0, ) —a;B; 


C, = fao fa(0) Cos(1, (0 -6,) ou C,= [de fa (0 +0,) Cos(1,0) 
S, -| dO f,(0) Sin(4,(8-8,) ou S,= fao fa (0 + 0,) Sin(2,6) 


(a,4,C, - B,S,) Fr 1 A 
T(r,0)= À,Cos(2,(0 -0,))- B,Sin(, (0-6, 
( ) È CAO fe | { L, | (a os( ( )) B in( ( )) 
La 
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Exemple : Secteur de disque annulaire (r,Ÿ) soumis à des conditions aux limites mixtes de Robin 
inhomogènes en r, homogènes en Ÿ 


Pour le problème général : 
1 à OT(r,0) 1 0° T(r,0) 

r + = 
r Or or r 0’ 
aT;(r,0)+ BT(r,0)., =0 
oT, (r,0)+ BT(r.0),, = fi(6) 


asTo(r,0)+ B:T(r,0)),, =0 


0 





aT, (r,0)+ BTO), =0 


D'après des calculs déjà réalisés sur les fonctions propres polaires dans la prise en compte des 
deux conditions aux limites : 


Fonctions propres polaires = ® (0) = a,4,Cos(4,(0 -0,))- B,Sin(2, (0 —0,)) 

À, valeurs propres solution de Sin(4,(0,-—0, a, aa 4 +p; B,)= À,Cos(1, (0, —-8 )XB,a., - Ba.) 
Fonctions radiales 

R(r)=Ar"+Br" R'(r)=å, (ar — Br) 

a+ BRCO), =a, (at, -BI > + (a + 81, )=0 

Aad, 7 + Bla” J= Blai, la -pl ) 

R(r) cle, ph + apla + Bada” y 

T(,0)= Ya, lo, -p> bA laaa Bl h> O 


Dans le a fonctions sinusoïdales, les normes se calculent ainsi : 
eO = 

T0 -0 Na,4,Cos(2,(0, -0,))- B.5in(2, (6. -0,)} +(a2,5in(2, (8, -0,))+ B.Cos(2, (8, -8))) 
ali — (a,A,Cos(1, (6, 8, ))- B,Sin(2, (0, - 6, ))\a;1,Sin(2, (8, —-8,))+ B,Cos(2, (8, -8,)))-a,B, 


n 
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La condition aux limites en r=l» donne la solution suivante : 
dyah pu) 

B,=(22,,""+81,+) 

R(r)= AT + B,r ” 

R'(r)=2,(4,r-B,r #) 

a R'(r)+ BR), =a,2,(4,1," "8,1, "")+8,(4,1."+8,1,*) 





r=l,; 


CRA) + B;R(Tr) = Apn (21 + Pi” }+ B,, (8,1, 7* eu a E 





r=l,2 


Con = 4m PAR T pis }- B,, (81. * = cs") 


[40 0), 0) 


> a Co: (0) = o (0) => a, — á 5 = B, 5 — aÀ,C, m 
Ši CaO OE CEOE Cn] 























B, = [d0 f,(0) ©, (0) =a,4,C, -BSS, 


n n 
0, 


C, = Í dO f,(0)Cos(4,(0-0,) ou C,= fao fa (0+0) Cos(4,0) 


0, 
0, 6-6, 

S, = | d0 f(O) Sina, (0-0) ou S,= [d0 f,(0+6,) Sin(2,0) 
0, 0 


iey E CE) 








n=1,0 | 
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Problèmes homogènes dans la dimension radiale pour un secteur de disque annulaire 


Revenons donc aux problème suivant: 
1 0{ OT(r,0) 1 O’T(r,0) 
r + 
r Or or r 00? 
GE) 0) 
or 


PEL ro) +B, T(r,0) 
Or 


= 0 





— BT(r,0) =0 





r=l; 


= 0 





r=l,2 


TEO =fr) 
T(r,0),.. = f,,(r) 


Et développons un peu plus le choix =A lors de la séparation, nous avons : 
r (r R'(r)+A Rr) =0 => R(r)= A'r’ +B', r” = A Cos(ALog(r))+ B„Sin(ALog(r)) 
©" (0)-2 00) = 0 => (0) = A', e + B'e% = A,Cosh(20)+ B,Sinh(10) 
L'équation différentielle radiale s'écrit : 
r (r R'(r))+uR(r)=r °R" (r)+rR'(r)+ u R(r) 
Posonsr = e' = dr = e'dt = rdt 
dR(r) __1dR() d°R(r) 1 dR( el P'RO) 
dr r dt dr’ r? dt r dt 
r R'(r)+rR'(r)+ u RO) =0> R" O+uRA=0 











R) = 4e" YBa E avec u> 0 ou u< 0 ou u=0 
Les solutions dépendent du discriminant de l'équation caractéristique de l'équation différentielle 


2 
linéaire soit 7” +4 =0>A=-—44 et s'annule à la valeur p=0. La valeur en p=0 donne : 
u=0=> R(t)= A4t+ B, => R(r)= A,Ln(r)+B, 


Exemple : secteur annulaire avec conditions aux limites de Dirichlet 


AT(r,0)=0 








T(r,0),., =0 
T(r,0), = fr) 


Les fonctions propres radiales sont : 


Si Ei 
P,(r)= sure Lol | Ze) = sites!) À, =nx La | 
x fl rl rl 


©” (0) = Sinh(À,0) 
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T(r,0)= $ 4,®(r)Sinh(2,0) 
Avec la solution se développant sous la forme d'une série : n=l 
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Les fonctions propres radiales sont orthogonales avec le poids w(r)=1/r : 


|= fu lormase-loies.. 
so E | 


D’(r)= Su Log Z 


rl 


En effet, avec les transformations de variables suivantes, il vient : 


La 
À. = rs Log Le > À, = ne. pe far z Sin, À,Log | Sin À,L0og T 
(le) li a pl r Li La 

Log F2. rl 


rl 








La Log(l,:/1,1) 
=> I = Í dr l Sin(à,Log(r"))Sin(4„Log(r'))= far Sin(2,t)Sin(2,,t) = 


1 F 0 


1 
=q [dt Sin(nr t)Sin(mx r)= 0 Le 
0 


n,m 2 

l 

Log| + 

mof- 
2 

On peut également se convaincre de l'orthogonalité en appliquant directement les résultats d'un 


système de Sturm-Liouville : 


r (r R'(r)'+u R(r)=0 => ar LRO) E Ror) =0 et u=X 
dr\ dr r 


1 
=> Système de Sturm — Liouville w(r)=—  p(r)=r  s(r)=0 
r 
Dans ces conditions, la solution du problème se calcule directement par la série suivante : 


T(r,t) = D A,D;(r)Sinh(2,0) 
n=l 





n 


[o Sink(2,8,) 


B, = far 20 sin À,Log ia 
r r I 


rl 





mens = y Be Sa à Logl | Saha O) 
F T | l Sinh(1,0,) l 


rl 
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On peut transformer les intégrales, afin de simplifier leur calcul : 








r= ser sr= / a = Lo( 2] 
La koef | La 
La 
f, (r ) Le f. 0 Log (1, 2/1 ) 
B, = far r Su Log -))- Í dr 222 Sin(A,Log(r'))= far Sin( 2t) f =a Jar Sin(nxt)f,.(r) 
ri K la 1 
avec À, = die 


I 
Log| > 
«| La | 


Soit une solution : 


TD > — BA one ISA ON a Lo 
koe! 2) &— Sinh(2,0,) be L 





r 


L, 1 





1 
B, =a C, <C, = | dr Sin(az t) fy, (2) 
0 


C, 
T(r,t)= Dre Los) jma, 0) 


n=l 


e koef | bodl 7 | 
es D . La La 
C, = [ar Sin(nx r)f,.(r) T 
0 


QAT 


= ©r=l e 


I a rl 
74 
La 


Lorsque la fonction limite est constante, il vient : 





í : 1 : 2 
fa (t)=1>C, 1e Sin(n x t )= (1-0) }= Cno Cra rT 


su auzoe[ 7 )\sm0,0) 
+00 C: , 4 T p 

Tr, =2Y Tr Sin 2,1 Sinh(2,0) = 

"o DAE a oef) HE Qn+1)Sinh(1,0,) 





n= 
nī 


D EN 
Log| => 
(ie) 


avec À, 
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Généralisation des conditions aux limites homogènes radiales 
Les fonctions propres s'écrivent sous une forme plus générale en revenant à la notation R(r): 


R(r) = À Cos 2 70) + UE 70) 


l t 
a= ol | r= os | E A relie dur r=l e" 
rl rl Lo! 2) 


rl 





R(t)= ACos(2at)+BSin(2at) 


Cela permet d'envisager d'autres cas de conditions aux limites homogènes comme par exemple 
une condition de Neumann de part et d'autre : 
dR(I,;) _ a dR(O) daR() _ 











Conditions ARa) =0 et 0 0 
dr dr dt dT 
Posons t= Lea - | t =Log(l.) t,=Log(l.,)—7T= 3 i aT= koef n | > dt = — 
rl 2 "1 rl 


R(r) = A Cos 2 70) +B sufa 70) 


R(r)= ACos(2at)+ B Sin(Aa r) 


R. (r) = sn = Eo 


— Conditions R: (t) = 0 





R: (t) = Aa (— ASin(A a t)+ BCos(Aart)) 
=> Conditions R.(0)=R (1) = A=B=0 sauf si B=0 et Sin(Aa)=0=>Aa=nr et R(r)= Cos(2a T) 


r 
L er 
(x) 
l 
L “r2 
oe Li | 


La norme se calcule de la même manière : 


| = far caf aroe E) = a [dr Cos(2,a r)f(T) avec À, = 
li k 0 


rl 





=> R,(r)= Cos| nm 


nī 


I 
Log| “2 
«| La | 





R, 








SO=R,()=Cos(2,ar) 





R 


n 


| 














: koef 12) 
i = a [dr Cos(,a rt) = 5 E 
0 
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Pour des conditions de Dirichlet en bas et de Neumann en haut : 





Conditions R(L,)=0 et RE 0 R(0)=0 et aRQE 0 
dr T 
R)= acol azos £ )) sso azoef Z) 
rl rl 
R(T) = ACos(Aat)+ B Sin(Aa r) 
R: (t) = Aa (— ASin(iat)+BCos(Aart)) 
=> Conditions R(0)=0 > 4 =0 
(2n +1) 


R(1)=0—=B=0 sauf siCos(ia)=0= 2a = m et R(T) = Sin(Aa r) 








La norme se calcule de cette façon : 


(2n+1)r 
l 

2Log| 2 

| La | 


2 2 «à! a Sin(1,at) a 
= a | dr Sin(à,ar) => [dr (Costa) = 5 T = = 
0 0 0 


R 


n 











1 
a far Sin(2,@ T)f(x) avec À, = 
0 


SO=R,(®)=Sin(2,at ) 


R 


n 





=>] 











À, 2 











Pour des conditions de Dirichlet en haut et de Neumann en bas : 
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Conditions Le =0 et R(l,)=0<& —=0 et R(0)=0 
r 


R(r)= ACof à Log 2 ) +B si 220812 )) 
r r 


Posons t= Log(r) t =Log(l.) t,=Log(l.,)-—T= GEN ys t-t = Loef a) 
t rl 


2 1 


dR(1) 
dt 





t, -t l 
>r=t in -ar ear = Log t2) 
a r 


R(r)= ACos(2ar)+BSin(1a rt) 

R(r)=2a(- ASin(iar)+BCos(2ar)) 

=> Conditions R(0)= 0 => 4=0 

(2n+ i) 





R@=0>B=0 sauf si Cos(Aa)=0= Àa = et R(r)=Sin(Aart) 











I r 
o) yale lE 
>R (r)=Si (2n+1)z cel r ) =( lc (2n+1)z Ly 


4 in A os 3 
Log la Log Le 
Li La 
La norme se calcule comme d'habitude : 





(2n+1)r 


2ro( 12 | 
La 


R 


n 











1 
a [ar Sin(2,a t)f(x) avec À, = 
0 


f(G)=R,(r)=Sin(i,ar ) 














P 1 
> |R, a farsn(aar) -4 fara- Cos(1,a tr) = SARS Sin(A,ar) _a 
2 2,a -2 
ft 
= |R,f == 
2 











Pour des conditions les plus générales de la forme : 
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Conditions di R, (La) ma B, RL) m 0 œ; R, (La) + B; R(,:) z 0 


rl rl 


R(r)= A Cos 2 70) + B si to )) < R(T)= À Cos(2a r)+ B Sin(2a T) 


Ro- 





S a Rl- ARU) =0 a, EO- g, RO =0 
(04 


rl 








S a, 0 D+ B: RG) © a, EO + p, RO =0 
(04 


R(r)=2a(- ASin(iar)+BCos(Aart)) 
= R(0)= À R(1)= ACos(2a)+ BSin(2a) 
= R(0)=Bia R(1)=4a(- ASin(2a)+ BCos( @)) 





g E E T EE E La 
al, La a À 
R Aa (— ASin(A a) + BCos(2a)) 


+ B,(ACos(2a)+ B Sin(2a ))= 0 








a, + B, RD =0> &, 
& r2 al, 


= à, À (-a,A Sin(2æ)+ Bl Cos(2a))+1,,B,(a,À Cos(2a)+ BL, Sin(2a))=0 
© Sin(a)ae,À - BBl,l,:)= 2 Cos(2aYaB.l,, +,Bl,;) 
R(t) =a A Cos(iart)+Bl,Sin(2at) 


Les fonctions ont donc la forme suivante : 


R, (r)= aaco atoe Z) + Pila so azos Z) 
rl rl 


À, tq Sina, alaa, -PBa )= 2,Cos(i,a Nabal, +a,Bl) 
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Les normes se calculent de la manière suivante : 


R 


n 














1 
za far (4,Cos(2,æ t)+ B,Sin(2,a t )) 
0 
1 
= a [ dr (4, Cos(à,a t} + B? Sin(A,a tÈ +24,B,Cos(3,a r)Sin(,a r) 
0 


1 
= Sfar (420 + Cos(21,at))+ B, (1— Cos(24,æ t))+24,B,„Sin(24,a r)) 
0 


n n 


je Sin(24,æ) já 1- Cos(21,a) 
24, 


Jar Sin( (24, QAT 
24,4 


1 
[ar Cos(24,a t 
0 


2 5 +B, )+ i (4, - B? )Sin(24,a)+24,B,(1- Cos(22,a)) 


R 


n 


=> 

















E (4 sa he ne) (4? - 8,7 \Cos(2,a )+ 24, B,Sin(2,o) 
Donnant le résultat définitif : 


R, = ar iri J PLA a2 pi Kosa) + 20B A Sinla,a) 


n 














2, 4 SinAaNaa A BB) 4,Cos(2,a Napl, +a.B,) 
On vérifie que l'on retrouve les normes déjà calculées : 
Cas(Da, =0 fB,=1 a, =0 B,=1 
































R 2 
4, 4 SnAa)=0 OC) 5 2 Log] || Pl _ a 
La La rl 2 
Cas (2)æ=1 B=0 aæ,=1 ß,=0 
R 2 
À, tq À, Sin(A,a)=0 AO 2, Log ai - r 
À, Li (,) 2 


Cas (3)a, =1 B,=0 a, =0 B,=I1 






































R 2 
A, tq À Cos(A.a)=0— Sin (4 a)=1 R, (1) = Cos| À, Log rA IES -= a 
À, La (2,) 2 
Cas (4)a =0 B=1 a, =1 ß,=0 
2 
R 
À, tq À,Cos(i,a)=0— Sin (4a)=1 a -sil 2 Loi] 4 J >75 
rl rl rl 
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Calcul des normes à l'aide de la théorie de Sturm-Liouville 
Avec le système de Sturm-Liouville : 


ERORO =0> (r r] ER =0 et u=% 
dr\ dr r 


1 

=> Système de Sturm — Liouville w(r)=— p(r)=r s(r)=0 
r 

Les normes des fonctions propres se calculent avec la formule : 


PD) DO à 4 reo] 





J dr w(r)® „(rọ = zef RE Ouôr 


ba 


On rappelle d'autre part que la valeur propre de l'opérateur de Sturm-Liouville dans notre 


__ 92 
problème était transformée comme suit : # = À 


En conséquence les dérivées paramétriques sur les valeurs propres changent de nature : 


u=% > du =2vdv 
Appliquées ces formules à notre problème, il vient : 


== pe)=r ®,()—R,(r) 





La 2 2 La 
pA R, (r) _ 1 K OR, (r) OR, (r) RO Ô R, (r) 
L, r 22 Or 0A OAOr 


ba 


D'autre part les fonctions ont la forme la plus générale suivante : 


R,(r) = A Cf à 70) + B Si (7) 


En normalisant l'expression des fonctions propres à l'aide de la variable adimensionnelle t, come 


ceci : 


r 





Posons t= Lea | t=Log(l.) 1, =Log(l,)—T= 


rl 2 1 


R,(r) = TE 70) + QUE 70) 


R,(r)= ACos(2ar)+ B Sin(2a rt) 
ƏR, (r) OR, (T) 
or 


, OR, (T) 
t R = 2 
RES ôl 


Posons R,, (r) = et R, 1(T)= 


RO RG e7 
ro ar a la 
R, (r) _ e @R, (T) 
0Aðr al, OÀ0T 
Ceci donne l'expression de la norme en fonction de la variable t : 


de | SR, | 
|R œ) ETI (roro R, (r) 2181 ji 


On a R,,(r)= Ri (©) 








De même 





aT= Lee 





r 


rl 


ja-% r=l e" 
r 
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Il ne reste plus qu'à exprimer les dérivées partielles, dans le cas le plus générale : 
R (t)= ACos(2ar)+ B Sin(a tr) 

R,.(c)=2a(-ASin(iar)+BCos(iat)) 
R,,(T)=at(-ASin(aar)+BCos(Aat)) 

OR,(r) _ OR. (r) _ 


0AT oÀ 
Après une petite session de calcul sur Mathematica, il vient : 


aa) 





a(-ASin(iar)+BCos(Aart)})}-2a°T(ACos(2at)+BSin(2at)) 


AB+2a(4° +B°)- 4 BCos(22 a)+ Sin(2à &) 





IRo = = 


_ AB+Aa(4°+B°)-AB(-25in(2 a) )+ (4? - 8° }sin(aæ)Cos(1a) 
2A 


>|R G = SA + B°)+ Sue) (4 - B° \Cos(1a)+24BSin(1x)) 


Ce qui donne le même résultat qu'avec le calcul direct des intégrales. 
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Exemple : secteur annulaire avec conditions aux limites de Dirichlet 





AT(r,0)=0 
T(r,0)\,.,, =0 
T(r,0),., =0 





T(r,0),., = f(r) 
T(r,0),, =0 


La solution se construit de manière similaire : 


reoat D Cu Su auzoe £ \sm(,(0-0,) 


x = Sinh(2,0,) 





rl 


C, = far Sin(nt)f, (T) 


kod 7 | 

La r 

T=T I = akoe =) 
koef | i 


rl 





La solution finale du problème : 





AT(r,0)=0 
T(r,0)|., =0 
T(r,0),.,, =0 





Tr, O) = f,(r) 
T(r,0), = fC) 


est donc : 


+00 


red? D Cin jS akoe £ ) |sat 0,)+Ž > Con Su auoe £) |si o) 


m Li Sinh(,0, P m Ki Sinh(À,05) 








rl 


C,, = far Sin(nt)f.(r) C,,= far Sin(nt)f, (T) 
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Continuons, avec une fonction limite linéaire : 


f()=r f,(r)=0 
T Loel 2 T Loel 22 
r=l,e" 4) > fre fe a) 


. La solution est : 





Sin 3 Log ) simn À,(0 -0 
+00 nl o L, (- 1)’ ] | La ( ( 0) 
= i Sinh(1,0 

| nr" + [Loge ) S 





Exemple : cas général secteur annulaire avec conditions aux limites homogènes mixtes de Robin 


Il s'agit du problème : 
OT(r,0) 106T(r,0) 1 O°T(r,0) 
+ + 
ôr? r o r? 00? 
CL. auT,(r,0)-BAT(r,0)., =0  a,T,(0,2)+ B:T(r,0),, =0 


ag T, (r,0)-B T0), =—f 0) AFT, (r,0)+ BETE, O) = fr) 
Problème qui se décompose en deux sous-problèmes : 
CL. a, T(r,0)- BAT(r,0) , =0 a;T(r,8,2)+ B,:T(r,0) , =0 


=0 T(r,0) fini Q=fr/1,<r<l,}x{0/0<0<06, <2x} 





at, (r,0)- BST(r,0),,=0 APTO + BPT EO) =R O 
et 
C.L. a, T.(r,0) —B,T(r,8). = 0 a, T.(r,0,z)+ B,:T(r,0) 0 
aëT.(,0)-BIT(r,0), =- 0) al ,0)+ B TE, O), =0 


Les fonctions propres radiales s'écrivent comme suit : 


a = Log l> at =Log| — 
La La 
R, (r) = Q hn cod aoe |] + ptasi atoe E) 
rl rl 


R, (t) = dd, Cos(2,a c)+ B,il,nSin(à,@ T) 
À, tq Sina ala 2, (2, ji _BaBlub)= 2,Cos(2,a Xa n B2l,, + a,>B,il,1) 


Avec les normes suivantes : 


|R, oF = Z (a, (A, y + Boeh 





Sin(2,a ) 


23 (an? (4, y m Bala XCos(4,a)+ 2a, BnA,l Sin(A ot ) 
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La transformation des intégrales deviennent : 


t 


r p É 
r'=— >e =r'>Tt= a = Log 
Li Log La La 
La 


R, r) = QA, cos aroe 7) + ZEZA) 
rl rl 


R,(r)=a,,À, Cos(A, t L Sin( Aar) 


PENI 











Log(( ralla) 


B= ju FOR Talon far O8, 0- a jar RORO 


ce qui permet de calculer les coefficients du développement en série et la solution complète du 
problème : 


a= Loef l> | AT -roe Z) 
La La 
R,(r)=@, À, cos aroe E) + BL. su[atos( 7) 
rl rl 


R,(r)=a,,À, Cos(,at)+Bl.Sin(A,at) 
À, tq Sin(2,e Ye, a, (,Ÿ = FERAS À ,Cos(4,a Xe iB,2l,2 + à Bali) 
Sin(À 
un Tla, (a „aY +B, 1 7) Ge), ‘(a aY -Bxl La 2 \Cos(2,a)+ 20,8, Al Sin(A, a) 
,,(0)= a; 2,Cosh(2,(8, -0))+ B; Sinh(2, (0, -0)) 
©, (0) =al, Cosh(4,0)+ Bê Sinh(2,0) 
4, = (2, Cosh(A,0 lat B? +a? BE )+ Sinh(2,0, Nate? (2,) + B2B?) 











Lo Laz 
Coordonnées r = C, „ = far R) fE) Cn = far RO ge (r) 
La r ln í 


1 1 
Coordonnées t = C, =a f dT RŒ) SYE) C, =a fdr ROS?) 
0 0 


Tr,0)= X 1 -(C,,8,,(8)+C,0,, ol a. co atos) + Bal, so aLes( £) 
F n n K rl rl 


Considérons maintenant la dégénérescence de ce même cas avec une valeur propre radiale nulle, 
r — 

soit Aar = de alors la solution est de la forme A+B Log(r), et la fonction angulaire de la forme A+B 9, 

coïncidant avec des conditions aux limites homogènes radiales de Neumann de part et d'autres. Ce 

qui conduit à une fonction radiale triviale : R(r)=1, dont on doit calculer la norme également avec 

le poids w(r)=1/r, car cette fonction propre est orthogonale avec les fonctions propres radiales . 


R,(r) 











, comme il suit : 
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Si di =, =] Ba = B,3 = 0 


Rr) =] R, (r) -cos aoe £ 


rl 


) R(r)=Cos(2,at) À, tq Sin(A,aæ)=0 


(R (r), R,(r)) = far RO R, (r) =a fdr R,(r)= a [dr Cos(2,ar)=0 
la r 0 0 


La solution pour ce cas particulier donne donc : 
Si di =, =] Ba = B,: =0 


a = Lo | AT -roe Z) AÀ, tq SinQ,a)=0— À, = “m 
a 


rl rl 





R,(r)= Ca” z Los! a J R,(tr)=Cos(n xt) 
a 


rl 


RO =S 


Oon (0) =afA,Cosh(à, (0, -0))+ B'Sinh(2,(8,-0)) ©,,(8)= a$2,Cosh(2,0)+ B? Sinh(2,0) 
Oo o(0)=a; +B; (0-0) 6,,(0)= a +230 
4, = (2,Cosh(2,0 at B? +a? Be )+ Sinh(,8,atat (à,) + Beee) 
4 = (a? B? + Bra? + B?B/0,) 
7... Cos TT Log _ ıı, Cos PT Log = 
$ a La 0 N Q La 0 
J Ji C) 
r r 


Coordonnées r = Ci, = [ar or) Cn far 
ba 


La 








2 Pdr L, 
R r)=1_| RO] aj e Pje 
la rl 








1 1 
Coordonnées tT => C,, = af dr Cos(n nt) @) Cp = af dr Cos(n zt)f (x) 
0 0 


Lz o Lz 0 
r r 
ae jar ee j'ar#-02 ) 
La La 





1 2" nm r 
T(,0)= (C0 (0)+ GB, (0)) +Z Z, (Cor Bon(0)+ GB (03) Ca "2 toef) 
0 nr=0 


n rl 
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Exemple: Représentation intégrale de la solution du problème de Dirichlet inhomogène sur la 
frontière d'un secteur de disque de longueur infinie. 


Bien que les fonctions radiales d'un secteur de disque plein ne soient pas apte à représenter une 
solution sous la forme de série sur un le secteur de disque plein, du fait de sa rapide oscillation 
autour de l'origine, on va s'inspirer des travaux de N Lebedev dans son ouvrage « SPECIAL 
FUNCTIONS AND THEIR APPLICATIONS, Prentice Hall 1965 », réalisé sur une configuration de coône 
sphérique. N Lebedev utilise ces fonctions qui ont une nature tout à fait similaire à celle du secteur 
de disque. On peut alors tenter d'exprimer sous forme intégrale la solution du problème de 
Dirichlet suivant : 

AT(r,0)=0 Q= {(r,0) € [0,+0]x[0,0, 

TO =0 Tr, O) o, = fa, (r) 

Condition 

Lim T(r, 0)= 0 Lim, fo, (r) = 0 


Comme nous l'avons vu les fonctions suivantes sont solutions de l'équation de Laplace : 
D'(r)= (a: Cos(r Log(r))+ B? Sin(r Log(r))) 

D°(0)= A? Sinh(t 0)+ B?Cosh(r 0) 

La condition en Ÿ=0, implique que cette solution se développe sous la forme : 

T(r,0)= (a: Cos(t Log(r))+ B} Sin(r Log(r)))Sinh(r 0) 


Supposons maintenant, comme le fait N.Lebedev, que la fonction limite sur la surface du cône 
admette un développement en intégrale de Fourier généralisée, soit : 


Posons Ë = Log(r)= dé = a 
7 


ROEI Í dr (F.(r)Cos(ré)+ F,(r)Sin(té)) 


© 


— = à fa Jo, (£)Cos(ré)= 1 fa (42 ) Cos(z Loc) 


= a dE fy (E)Sin(ré)= 1 fa [40 a Cr sule Lost) 
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Le respect de la condition aux limites de Dirichlet conduit à la solution suivante : 
T(r,0)= (a: Cos(t Log(r))+ B} Sin(r Log(r))}Sinh(r 0) 


Développement intégrale sur le paramètre t 

=> T(r,0)= Jala: Cos(t Log(r))+ B? Sin(r Log(r))}Sinh(t 0) 

Comme TEO, = fo, (r) 

far(a! Cos(r Log(r))+ B} Sin(r Log(r (r))}Sina(r 0 )= Jar (F, t )Cos(té)+ F.(r)Sin(ré)) 
0 


go Ee p RO 


| Sinh(r 0) ‘  Sinh(t 0,) 
T(r,0) = arr. (e)Cos(r Log(r))+ F(t )Sin(t Log(r pse o) 


avec 
= ES jo (22 Cos(r COR = -1 fa a [0 2 Sin( T Czego) 
m 0 4 0 


Dans la littérature, il y a une condition pour que l'intégrale de Fourier existe sur la fonction limite 
f(r) (voir E. C. Titchmarsh, Introduction to the Theory of Fourier Integrals, second edition, Oxford 
University Press, London (1950), Théoreme 3, p. 13): 

Jo, (r) continue 


Jo, (r) bornée sur tout interval [r1,r2] 


far Li FAQ) 


Avec le problème complémentaire 
AT(r,0)=0 Q={(r,0)e[0,+x]x[0,6, |} 


TEO ERO TEO =0 
Condition 
Lim, T(r,0)=0 Lim, pofo r)=0 


La représentation intégrale de la solution est : 


T(r,0)= Pre (r)Cos(r Log(r))+ Fè (t)Sin(t Log(r))) ee) 





est finie 





avec FX = je (£2 œ RO dosi Lot) F(r)= 1 far (£2 Sin(r Lost) 
mS r 
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ce qui fait que par principe de superposition pour le problème de Dirichlet inhomogène sur la 
surface du cône : 


AT(r,0)=0  Q={(r,0)e[0,+]x[0,0, |} 

T(r,0),, ENO TEO o, = Ja) 

Condition 

Lim,,,.T(r,0)=0 Lim ofo r)=0 Lim, fs (r)=0 


La représentation intégrale de la solution est donc la suivante : 


T(r,0)= far; (t)Cos(t Log(r))+ Fè (t)Sin(t ie (0, 0), 





Sinh(r 0,) 
+ Jar E Coste Loue) Fie) Sin Los) SE D) 
avec 
Fi(r)= 2 Es At Cos(r Lot) Fi(r)= 2 < (£0 Sin(r Lot) 
m r m r 
A = fo (22 JO cos( EL) f; =— -ifa (42 JO Sinf Czesc) 


avec les conditions suivantes 

Jo (Y), fo, (r) continue 

Jo Y), fo, (r) bornée sur tout interval [r1,r2] 

TO Tpha 

[a ; | et [ar "i | 
r 


0 0 





sont finies 
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Extension au cas cartésien avec l'angle d'ouverture du cône à x 


Lorsque l'angle d'ouverture s= n, le problème se développe par exemple dans le demi plan 
supérieur : 

AT(x,y)=0 Q = Kx, yje [- %,+00]x [O, +00 |} 

ÉD EE AO C0) PET AC) 

Condition Lim, _..,,T(x, y)=0 Lim, „fè (x) =0 Lim, „f> (x)=0 


J—+o 


La représentation intégrale de la solution devient 
Si. {x,y}e [0,+c0[x [0,+c] 





D), 





+0 Sin (= — Aretan 
T(x, y) = far(r:(e)Cose Logl Je +y )+ F: (r)Sinle Log =y”) 


Sinh(r x) 
Y 


3) 


Sin Aretan 


QE Areta 











+ fadle: (r)Cosle role +y’ )- FE (c) Sinke Logli j +y )) 


0 Sinh(r x) 
Si {x, y} € [- %,0]x [0,+] 





X 
X 








), 
Sinh(r x) 


D) 


T(x, y) = far: (r)Cosle ZAN +y’ )- F; (c) Sink- Log(x7 +y”) 








, Sin Fr — Arctan 
+ fale Gene Lee rs ror en 
Avec r= +» 


0 = Aretan | 
(e E + Eak] 


{x,y}e }e[0,+o]x [0,+c] 
x 
F° (t) = 1° fa (a6 2 Cos(r Lot) F? (t) = e (A2 Sin(r Lot) 








0 = x — Arctan 











T 


Fr)-1 ja [re 


a ta 
a 





T X 


Cou Log()] F = fa a [Fe Œ E Sinf (Log) 


0 
et les conditions suivantes 


Jo Œ), fo, (X) continue 

Jfa (x), fa, (x) bornée sur tout interval [x1, x2] 

+00 r(x je 7 (x 

fa À. fax ) 
x 


0 0 





sont finies 
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Remarque sur le respect de la condition sur les fonctions limites : on a vu qu'il peut exister un 
développement intégral de la solution du problème aux limites si la condition suivante, est 
respectée : 

f.(r) continue 


f.(r) bornée sur tout interval [r1,r2] 


ja dr' A 





est finie 


prenons par exemple un problème d'électrostatique en deux dimensions consistant en une charge 
placée sur la bissectrice du cône à la distance a de l'origine au dessus de l'origine. Il s'agit d'un 
problème de Poisson où le terme source est constitué par la charge électrique. On peut vite voir 
que le problème se résout par la superposition de deux problèmes, l'un du potentiel d'une charge 
dans l'espace placé au même endroit, sans condition de limite fixé, et l'autre le problème induit 
par la conditions aux limites de valeur opposée au potentiel sur la surface du cône. Sans entrer 
dans le calcul exacte, on obtient un potentiel de la forme 


V(r) =- 2 Log +a? —-2ar Co 0 — 2) 


Ce qui conduit à des fonctions limites sur la limites du cône de la forme 
Jo = fa r) G +a°-—2a rcos( S )) 


Et bien cette fonction ne respecte pas la condition imposée, car le problème se pose tant dans 
l'intervalle inférieur d'intégration que sur l'intervalle supérieur d'intégration. Si nous prenons une 
borne inférieure C2 pour l'intervalle d'intégration supérieur suffisamment grande, alors : 


A E O A0). a) 
far = = [ar A : er 





r 


0 0 ‘4 


Voyons ja LA) ) prenons C, tq koer +a’ — 2ar Co 2) >1 V0, € [0.27] 





0 
T ILog(r? +a?’ —2a r) Les Lol HR caf) Fe Log(r° +a’ +2a r) 


far 7 < far 7 < far x 


C C2 C2 


koer +a°—2ar ca( à) 


r 




















Or far toeta | ge il suffit que C, >a +e' > far ni: ar 
AE r r 


à) Ca Ca 








Dans ces conditions le problème posée d'électrostatique en deux dimensions ne trouve donc pas 
de solution en représentation intégrale. 
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Prenons maintenant des conditions aux limites de la forme : 


y P r? sur[0,7] 
Jo Ee= far) °” B>0 
0 sur[n,+] 
nÍ U 
Commençons par le cas simple B=0, alors la condition devient: ° i , ce qui indique 


clairement que l'on ne peut envisager de représentation intégrale de la solution, quand bien même 
dans ce document nous avons donné un développement en série de cette la solution avec B=0 
pour les mêmes conditions aux limites de Dirichlet. 
Prenons maintenant le cas B quelconque, alors la condition est remplie : 
To r. B 
R] drr?” = finiesiß>0 
i p | 
On va donc calculer la représentation intégrale : 
RO=ROSK()=E()=F() e F()=Fi(r)= FT) 


T(r,0) = fa (F“&)Cosl ere RE Log(r) (Sinh(r (0-0, ))+ Sinh(r 0) 





avec 


F(r)= Tig rP Cost Log(r)) Fe(r)= Tig rf Sin(t Log(r)) 
m 0 m 0 
To Šo Šo 
far rf™Cos(t Log(r)) == fag e” Cos(t €) = Les of. -p fag e” Sin(t :) 
0 —00 T —00 


ro ĉo So 
far rf" Sin(t Log(r))= S= fag e” Sin(t £)= z- le” Cos(r EE + B fag e” Cos(t 5) 
0 —00 T =% 


E E E 


T 
sfr: E) Š (- le” Cos(r &)Ë, + B fe sin(e of. ) > S- Erle cos(e DL. = ee 33 
T T T G +8 
(8 e Sin(t &,)-re/* Cos(r £) E C= est EJop (B e° Sin(t é, )- re" Cost &) ) 
(+p?) T x t+?) 
(8 ee Cos(r é )+e®®r Sin(r £, ) 
(z? + B°) 
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En insérant ces expressions dans la formule intégrale on obtient la représentation intégrale de la 
solution du problème aux limites : 


(8 e”™ Cos(t é )+ er Sin(r £, )) s(r)= 1 (8 e” Sin(r é, )- te”™ Cos(t £) ) 
m (+p?) F ( ) z (+p?) 
(F°(t)Cos(t Log(r))+ F* (t)Sin(t Log(r))) 

— eP5o ((8 Cos(r ë)+t Sin(r Gi ))Cos(t Log(r))+ DE E)-T Cos(r Čj ))Sin(t Log(r))) 








F(c)= 





m t +B’ 


ef (B Cosle Log(r)-rE,)+r Sin(s Logfr)-58,) _r G Ga: 2er) sn 26) 





mT (t? + B°) m (t? + B°) 


se L cof: Leaf} Sul tr) (Sinh(r (00, ))+ Sinh(e 0) 


T(r,0)= a] dr Cp) Sinh(r 0,) 





Pour B=1, soir f(r)=r, nous avons : 


i (cof: 2H) Sul 74) (Sinh(e (0 -0,))+ Sinhle 0) 


2 To 
SYS T 14e (r? + 1) Sinh(r 0,) 





Prenons maintenant des conditions aux limites inspirée de potentiel électrostatique induit par une 
charge ponctuelle, placée à la bissectrice de l'angle à une distance a de l'origine. La fonction est de 
la forme : 




















(= L = NES = 
pta (2) y +a° — 2a rCos (2) 

PAOE T e daed a 
Z+- 2Cos (2) r?’ +a° -— 2a rCos (2) 
Ņ\r 2 \ 2 
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Le calcul de la solution intégrale passe par celui des intégrales suivantes 








Fe(r)= : Í 4 To Cos(r Log(r)) 
Tr [£ r (2) 
+ 2Cos 
r a 2 
a r dr 1 
En posant é = Log(r)=>—+ — = 2Cosh(£ — Log(a)) d£=— s=£-Log(a) de plus fas 
r a r 


7 Jeo- Co“ 
2 





Sin(r s) 
































0-2 fæ Cos(r é) - 2 nf Cos(t AEO T,Cos( -U ja 
NT Jeong — Log(a))- cf?) Fe eoo- co?) Co: 
F; (t = far [40 (r Le O Sile Lt) - asne ces) [as Cos(r s) 
me ki coss )- Cos (2) 
Fo, e fo [22 RO cos(r Lost] Te Crea fas CE s) 


de a Cosh(s)- co?) 


Fe 1 l Far (22 aO Sinf 20) - ae ogia) fa cts) 
o Cosh(s)- co ® 
2 








Les intégrales sont connues : 





















































[as Cos(r s) _ 2 [as Cos(r s) 6 = a Pa [- co( 2) 
5 Cesi(s)-cos( 2) i cosmo) co à) one 
F«(r)= T,Cos(r a) ja Cos(r s) = T,Cos(t Log(a)) „ -cos 2) 
i 9 Cosh(xr) Se 2 
Cosh(s)- Co à) 
\ 2 
Fi (t)= T,Sin(t cote) ja Cos(r s) 2 T,Sin(t Log(a)) , -cos 2) 
i 9 Cosh(xr) tir 2 
Cosh(s)- co?) 
\ 2 
3 T,Cos(r Log(la)) Cos(r s) _ T,Cos(t Log(a)) | 2 ) 
F, P — Cos| — 
OUR go de) OU 
2 
P T Sin(t Log(a)) Cos(t s) _ T,Sin(t Log(a)) E (2) 
RO 92 Cosh) l (2 


zR ET Cos (2) 
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Et la solution intégrale s'écrit finalement : 


P [- Co$ )) 
T(r,0)=T, [ar wir Cos(r Log(a))Cos(r Log(r)) 
»0)= 15 A Cosh(xt )Sinh(t 6,) + Sin(r Log(a))Sin(r Log(r) 





| sm (e,-0)+ 





+0 Fa en 2 os\t Log(a))Cos(r Lo + 
th [az rueda ent rasent sast y j» M 2 
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= T(r,0)=T, Í dr Co Loft )) Cda |- co( 2) 


Cosh(rr )Sinh(t 0 À ) a T 
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